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I . INTRODUCTION 


The majority of computer aided circuit analysis programs belong to the 

class of "numerical programs", that is, the output is some numerical value. 

(8)* (9) 

At the time of our research a few programs, most notably ANP1 V 1 , NASAP V , 

and CORNAP^ 8 ^ could generate network functions as rational functions of s 

but did not allow the value of any element to be left in symbol form. The 

research project reported here represents, we believe, the first effort to 

generate symbolic network functions. By a symbolic network function we mean 
V _ V _ 1 _ I 

out out out out _ . . r . . . . , , 

, - , — , or as a ratio of two polynomials of one of the 

in in in in 

following types: 

(1) all network element values are represented by symbols (the symbols need 
not all be different) 



Examples: 

V out 

s 2 LRC 


V in 

s 2 2LRC + s(L+R 2 C) + R 



V out 

ZYR 



V in 

2ZYR + Z + R 2 Y + R 

(2) 

some element values are specified numerically, some 


Example: 

V 

out 

s R 


V. 

in 

s 2 2R + s(.5x 10 6 +150R 2 ) + .75xl0 8 R 


or 



(3) 

all element values are given numerically. 


Example : 

V out 

2 

s 


V in 

2s 2 + 2 x10 4 s + .75xl0 8 


There are many reasons why one may be interested in totally or partially 
symbolic network functions. The following presents a few of the more important 


ones. 

(1) Insight. To illustrate the added "insight" symbolic programs can provide 
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in comparison to numerical type programs, suppose we have been asked to 
verify that the network in Fig. 1 is a negative impedance converter for 
large P, i.e., 

Z, (s) -» -Z T (s) as P ■* ® (1) 
in L 


X 



Figure 1 

To verify (1) with some degree of certainty using a numerical program 

would require evaluating Z (s) (and Z (s)) for many different values of 

in l 

P and frequency u), a time consuming process at best since most programs 
must completely re-evaluate the network response for every relatively 
large change in parameter values. Furthermore, the resulting verification 
would only be valid for the particular structure and component values 
chosen for Z (s). With a symbolic program, one computer run gives the 
symbolic transfer function 

Z i„ (s > ’ Z L (s) < H~S ) 

from which (1) follows immediately. 

(2) Error Control, To demonstrate how a symbolic program can be used to 

effectively control round-off error, consider the differential araplifier^^ 
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shown in Figures 1 and 2 of Section II-l. If network branches 1, 2, 5, 


and 11 are chosen as the tree for deriving the signal-flow graph (SFG), 
then the set of nontouching loops (all orders with sign attached) 
which belong to the numerator of the low frequency transfer function 



s=0 


N 

A 


is given by 

[R1R3A2 -R1R3A1 R1R3A1 -R1A1 -R3R1A2 R1A2| 

(RE) 2 R2 * (RE) 2 R2 ’ (RE) 2 R2 RER2 ’ (RE) 2 R2 ’ RER2^ 


and the corresponding set for the denominator is given by 

jRl R1 R3 R3 R1R3 R1R3 R1R3 R1R3 1 
lR2’ R2 * RE’ RE’ R2RE ’ R2RE ’ R2RE ’ R2&E / 

Letting A2=A1, R1=5K, R2=15K, R3=10K, and RE=25, evaluate N and A by 
summing the terms in the order given in the above sets keeping each 
number generated to 8 significant digits. Then 

N - A1 [5. 3333333-5. 3333333+5. 3333333-. 013333333-5. 3333333+0. 13333333] 
= 3.3x10 _8 A1 


and 

Thus 


A = 

I „ 
out 

V. 

in 


1335 


s=0 


3.3x10 

1335 


A1 


which is incorrect since N^. Although the above transfer function was 
derived using SFG theory, round-off errors which cause erroneous results 
can occur in any computer program restricted to numerical evaluation, 
and are generally very difficult to predict or control. Because round-off 


error enhancement in the evaluation of network functions often occurs 


as a result of widely separated values of some of the network elements, 
one method of error control would be to leave such element values in 
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symbolic form. This technique can be applied to the above example by 
noting that RE should be kept as a symbol since its value is considerably 
less than the other resistance values. Thus, keeping RE as a symbol and 
re-evaluating N gives 


3333.3333 _ 3333.3333 3333.3333 _ .33333333 

(RE) 2 (RE) 2 (RE) 2 ^ 


3333.3333 

(RE) 2 


+ .33333333 1 
* RE J 


0 


That is, 


I 

V 


out 

in 


s=0 


0 


(3) Sensitivity Analysis. Sensitivity of the network function to changes in 
a particular network parameter can be found using a symbolic program by 
keeping this parameter as a symbol and then performing the required 
differentiation. Although there exist powerful numerical techniques for 
sensitivity analysis, the above procedure using a symbolic program has 
the particularly desirable feature of being less susceptible to round-off 
errors. 

(4) Parameter Variation. Suppose we wanted to evaluate the network function 
for many different values of one or more network parameters. Using a 
symbolic program, we could leave these parameters in symbol fora and 
then efficiently and accurately perform the large number of required 
evaluations on the resulting symbolic network function. On the other 
hand numerical programs now available must re-derive the transfer function 
for every relatively large parameter change. 

(5) Iterative Piecewise Linear Analysis of Resistive Nonlinear Networks. 

Part of this powerful analysis technique requires the solution of a 
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resistive linear network where some resistances and some d-c sources 
are kept in symbol form. 

The primary objective of this project has been the development of some 
new or improved concepts needed to make a symbolic network analysis program 
efficient with respect to program storage and execution time. The project 
culminated in the program SNAP (Symbolic Network Analysis Program) which 
finds symbolic network functions for networks containing R, L, and C type 
elements and all four types of controlled sources. SNAP contains the following 
unique features: 

(1) The extensive use of a path-finding algorithm in place of matrix operations, 

(2) Efficient techniques for finding all loops of the SFG and for enumerating 
all higher order loops, 

(3) The use of the "compact signal-flow graph" instead of the "primitive 
signal-flow graph", and 

(4) A simple coding technique which is used 

(a) manipulate symbols thereby allowing the complete program to be 
written in Fortran (another important aspect of the coding scheme 

is that it permits repeated symbols to be treated as one symbol), and 

(b) determine whether or not loop sets touch in the algorithm for 
enumerating higher order loops. 

New techniques for handling multi-inputs and multi-outputs are also presented 
in this report although they have not yet been incorporated into the program 
SNAP. 

* (3) 

At about the same time the results of this project were disseminated, ' another 

symbolic program^) (by coincidence also called SNAP) whose primary concern is 

the on-line use for design purposes made its appearance. 
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II. A GENERAL DISCUSSION OF THE BASIC ALGORITHMS 

II-l . Fonsulating the Signal-Flow Graph (SFG) 

a. Data Required 

A SFG is generated by SNAP (Syabolic Network Analysis Program) from data 
specifying the topological structure of the network, the input-output variables, 
and the characteristics of each network branch. The network topology is 
described by 

(a) A unique number for each branch, and 

(b) the initial and terminal nodes of each branch as determined by the 
assigned current direction.. 

The input to the network must be a single independent voltage or current source 

and the output requested must be the voltage or current associated with a 

•k 

network branch or the voltage between any two nodes of the network. Finally, 
each branch is characterized by 

(a) a symbol which specifies its type, i.e., 

passive branches: R, G, L, C, Y, Z 

control sources: VV, VC, CV, CC 

independent sources: E, I 

(b) a symbol representing the branch name together with the branch value 
if specified, and 

(c) the branch number of the control (for dependent sources only) 

As an example, consider the network given below (from a paper by A. DeMari^^). 


Refer to Appendix A of section I1I-1 for a discussion on how to handle 
multi-inputs and multi-outputs. 
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Figure 1 Figure 2 

Table 1 (Network Data) 


Branch 

Type 

Branch 

Number 

Initial 

Node 

Terminal 

Node 

Symbol 


Value 

Control 

E 

1 

2 

1 



•» 


R 

2 

1 

3 

R1 

« 

5x10 


C 

3 

1 

3 

C 




R 

5 

1 

4 

R1 

■ 

5xl0 J 


C 

4 

1 

4 

C 


■» 


R 

6 

3 

4 

R2 

n 

15x10 


CG 

7 

3 

2 

A1 



10 

CC 

8 

4 

2 

A2 



9 

R 

9 

2 

5 

RE 

m 

25 


R 

10 

2 

5 

RE 

EB 

25 3 


R 

11 

5 

1 

R3 

“ 

10x10 



b. Finding a Tree 


The formulation of a SFG starts with the choice of a network tree. The 
selection of network branches to be used in the tree is made as follows: 
Independent voltage sources and controlled voltage sources are the first ones 
to be used. Then come the passive RLC elements in any order. In choosing the 
(J+l)th branch, the undirected graph formed by the J branches already selected 
is tested to determine whether a path exists between the two terminal nodes 
of the (J+l)th branch. If so, the branch under consideration is disqualified. 
If not, the (J+l)th branch is added to the tree. Let n be the number of nodes 
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of the network graph. When n-1 branches are successfully chosen by the above 
process, we have obtained a tree. 

As an example, consider the network of Fig. 2. If, following the selection 
of the voltage source, the passive branches are examined in the order by which 
they are listed, the tree shown in Fig. 3 results 

Tree branches: 

b l “ 1 * b 2 “ 2 * b 3 " 5 » b 4 “ 9 
Links: 

\ =■ 3 » * 2 “ 4 » l 3 “ 6 » \ “ 7 
«* 8 , “ 10 , t » 11 


Figure 3 

It is important to note that the complexity of the SFG and consequently 

the time required to evaluate the transfer function depends on the tree 

( 2 ) 

selected. A brief summary of the rules for choosing a "good" tree is 
given in Appendix C at the end of Section III-l. 

c. Rules fog Formulating the Compact SFG 
A "compact SFG" is a signal- flow graph whose node variables consist only 
of tree branch voltages and link currents except when additional nodes are 
needed for control sources or for the output variable. This type of SFG 
can be more efficiently evaluated than the so-called primitive SFG which 
contains one node for the branch voltage and another for the branch current. 

The compact SFG is constructed according to the following rules: (An 

example as derived from Fig. 2, Fig. 3, and Table 1 is given in Fig. 4). 

Rule (1): For each link the unique fundamental circuit containing 




Figure 4: SFG 


Initial 
Node 


6 

2 

3 
5 

4 


Table 2 (SFG DATA) 


Terminal 
Node 


Exponent 


Branch 


1 

3 
2 

4 

5 


Branch 
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branches bj,, i°l,2,...,m is found. Tiro sets of SFG branches can then be created. 

Set (a): For each passive branch in the tree branch set b^, i*=>l,2, . . . ,m, 

a directed branch in the SFG is fonsed from node I . to node V with 

b i 

weight equal to the impedance of branch b^, prefixed with the proper 
sign (positive, if the directions of ^ and b^ concur in C^, and negative 
otherwise) . 

Set (b): If the link is a passive branch, a directed branch In the 

SFG is formed from each node V, , i“l,2, . . . ,m, to node I. , having weight 

b i \ 

equal to the admittance of link 1^, prefixed with the proper sign 
(negative, if the directions of and b^ concur in C^, positive other- 
wise) . 

Rule (2): If any of the four types of controlled sources are present, a 

directed branch is created in the SFG from the controlling variable to the 

controlled source, having weight equal to the constant of proportionality 

(g , beta, etc.). If the controlling variable is a link voltage or a tree 
m 

branch current, one more node is added to the SFG to represent this controlling 


variable X (node Fig* 4 is a node of this type). X is then expressed 

in terms of the tree branch voltage or link current through a simple immittance 
relationship. 

Rule (3): If the desired output Y is neither a tree branch voltage nor a 

link current, then one node is added to the SFG to represent Y. Y is then 
expressed in terms of tree branch voltages or in terms of a link current 
through a simple immittance relationship. 

Rule (4): Finally, the SFG is "closed" by adding a branch with a symbolic 

weight (FB), directed from the output to the input node. 
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d. The Gala Formula for "Closed" SFG ^ 4 * 

The purpose for introducing the closed SFG is because only all orders 

of nonrtouching loops need be found as opposed to the evaluation of Mason's 

formula which required enumerating certain paths as well as loops. 

To derive the gain expression for the closed SFG consider first Mason's 

equation for the transfer function. 

m 


x 

_ O l®*! 

T o — — » ■ "■ ■ ■■■* 

x i * 


where 

A = 1 + E(-1)^EL, is the determinant of the SFG 

J th 

L. . i s the product of the transmittances of the k set of non-intersecting 

k, J 

loops of order j. 

P^ is the transmittance product of the i^ path between X^ and X^ 

A^ is the partial determinant obtained from A after removal of all 

tVl 

loops intersecting the i path between and X q . 

Let A be the determinant of the closed SFG. It is then noted that since 
c 

is the set of all paths from X^ to X q , the loops present in the 

closed SFG not present in the original SFG will be precisely C(FB)P^,}° a ^ 
where FB is the symbol assigned to the added branch. Further, since the 
path FB contains only nodes X^ and X q which, in turn, are present in every 
path F^, i a l,2,...,m, it follows that the non-intersecting loop combinations 

that do not touch the loops (FB)P^ , i=l,2 m will be precisely those 

combinations which do not touch the path P^, i“l,2,...,m. It follows that 

m 

A - (FB) E P A + A 
C i«=l 1 1 

Thus, the transfer function can be found by simply sorting the terms of the 
determinant of the closed SFG. 
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II-2. Manipulating SFG Branch Weights 

Each branch weight in the SFG is of the form 
Constant • Symbol ♦ 8° 

If an arbitrary branch has an initial node X^ and a final node X^, then the 
three parameters 

C(X^,X^) ® constant 
S(X i# X f ) =» symbol 

and E(X^,X^) « exponent of s 

completely define the weight of the branch. After a loop or a set of nontouching 
loops has been found in the SFG, say by some path-finding technique, it is 
desirable to combine the weight parameters of each branch in the loop set 
to fora a composite loop set weight. The loop set constant may be easily 
formed by taking the product of the constants associated with each branch. 
Similarly, the loop set exponent parameter is readily found by summing the 
exponents assigned to each branch. However, because computers are not 
particularly adept at symbol manipulation, it is inefficient with respect to 
both time and storage to form directly a composite loop set symbol. A much 
better technique is to convert each branch symbol into a numeric code. These 
codes are assigned to the SFG branches as follows: Each distinct symbol in 

the SFG is stored in the array S(j) and assigned a code B^ where B is some 
base Be{2,4, . . . , 2 C1 3. How for an arbitrary SFG branch having initial node 
and final node X^ which contains the symbol S(n), the code 

K(X 1 ,X f ) « B n 


is assigned. 
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The real value of this coding technique sterna from the fact that the compos- 
ite loop 8et code formed by summing the codes representing the Individual 
branch symbols can be uniquely decoded provided the number of identical 
symbols combined into any code is less than B. 

As an example of the above concepts for manipulating the SFG branch 
weights, refer to the SFG shown in Fig. 4. Consider, in particular, the 
loop defined by the node sequences 

v 2 - I 3 - V 2 and V 4 - I 5 - V 4 


Then 

3 3 

composite loop set constant ■ (-5x10 )(l)(-5xl0 )(1) 

- 25xl0 6 

and 

composite loop set s power a 0+l-S-0 + l“2 

To find the loop set code, an array of distinct symbols of the SFG and their 
corresponding codes must be set up. 

no symbol < ■> 0 

S(l) - FB < >4° 

S(2) = C < >4 1 

S(3) *= A1 < > 4 2 

S(4) “A2 < > 4 3 

Hote that because there will be at most two identical symbols in any code, the 
base 4 was chosen. Using the above codes gives 
composite set code » K(V 2 ,I.j) + K(I g»V 2 > 

+ k(v 4 ,i 5 ) + k(i 5 ,v 4 ) 

■ 4* + 0 + 0 + 4* 

0 8 


How to decode this number, say in the output, we would write 
8 “ 2(4*) 

- (S(2)) 2 



which is indeed the symbol associated with the loop isnaittance product. The 
above coding schema for manipulating symbols is easily adapted to the computer 
by incorporating the masking operation .AMD. . To determine the number of 
S(l) type symbols contained in a given code, the .AMD. operation is applied 
to the code and B-l. In general, the number of S(J) symbols is found by 
dividing (using integer division so as to truncate the remainder) the code 
used to determine the number of S(J~1) symbols by B ©nd then applying the 
.AMD. operation. For example, consider the loop set previously discussed, 
loop set code “ 8 “ (000000000100 ) ^ 

B -1 - 3 » (000000000011) 2 
(loop set code) .AMD. (B-l) “ 8. AMD. 3 

“ 0 

Thus, the symbol S(l) *= FB is not present. Mow divide the loop code by B 

and repeat the above procedure 
@ 

new code <=» -r ® 2 
A 

(new code) .AMD. (B-l) “ 2. AMD. 3 

« ( 000000000010) 2 
» 2 

This Implies C is contained in the code 8 and that its exponent is 2, i.e. 

2 

C . The process stops when the code is reduced to zero. 

Each loop set (of any order) contributes to a terra in the network function. 
As each loop set (of any order) is generated and coded, it is compared with 
existing terms. If a term with the same symbol code and power of s exists, 
then the constant of the term is updated by adding to it the constant of the 
new loop set. Otherwise, a new term is created. Note that the above process 
is an Important step towards reducing the storage requirements. 
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After all loop secs have been found, the transfer function Is complete , 
and it retrains only Co transform Che symbol code of each term into its corres- 
ponding symbol set by the .AND. operation previously described. 
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II-3. Generating First Order Loops 

a. General Description 

Let the nodes of the SFG be labelled 1,2,..., N. All first order loops which 
contain node J (J®1 Initially) can be found by conceptually splitting node J into two 
nodes, one node containing all incoming branch and the other containing all 
outgoing branches, and then enumerating all paths between these two nodes. 

All branches going into node J are then removed and the process repeated for 
node J+l. Clearly, this procedure will produce all circuits with no dupli- 
cations. 

The problem of efficiently finding all circuits now becomes one of finding 
paths. The path-finding algorithm utilized by SNAP is based on a routing 
technique which conceptually resembles that proposed by Kroft^^. However, 
because our ultimate objective is a flexible user-oriented program, we have 
chosen to use FORTPJVN instead of SNOBOL as Kroft did. A general description 
of the concepts contained in the algorithm will be given here in addition to 
a rigorous step-by-step description presented at the end of this section. 

Consider the SFG given in Fig. 4. The topological structure of the SFG 
can be completely described by the following routing table where the entries 
in the row are the set of all nodes of distance one from node J. Note 
that the entries of each row are made to decrease as the column subscript 
M increases. This facilitates modifying the table after all paths through 
a particular node, say node J, have been found because only the right most 
non-zero entry of each row must be tested, i.e. if that entry equals J, it 
i 8 a set to zero. As an example in using the routing table, the following 
two circuits can easily be shown to compose the complete set of circuits 
containing node 1. 

1-11-9-12-8-5-6-1 


1-11 -9-10-7-2-6-1 
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R(J,M) - 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


11 

6 

2 

5 

6 
5 
2 
5 

12 

9 

9 

8 


3 


4 

2 1 


11 10 
7 


Routing Table 

A particularly important feature of the path- finding algorithm is the 
method by which each new node generated from the routing table must be tested 
to prevent loops from being formed. Rather than comparing the prospective 
node to each node already in the path, it is much more efficient to define 
the function 


F(X) - {J 


1 if X is contained in the path node sequence 
if I is not contained in the path node sequence 


on which only one logic test need be made. 

Additional insight may be obtained by viewing the path-finding technique 
graphically. That is, the process by which paths are generated can be 
observed by applying the following two rules directly to the SFG. 


(1) Let node J be the last node added to the path node sequence (initially 
J ■ input node). To select the next node, traverse that branch 
connected to node J that goes to the highest numbered node satisfying 
both the following requirements: 



1U 
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a. we did not Just back up from this node while applying rule 2, 
and 

b. this node is not included in the path noae sequence. 

Repeat this process until the output node is reached (then store 
the node sequence and go to rule 2) or until no new node can be 
found the satisfies (a) and (b) (then go to rule 2) 

(2) Back up along the path just found (this is always possible unless 
we are at the input node in which case all paths have been found) 
until a new route can be taken according to rule 1. 



Hie heavy lines of Fig. 5 show the path which results from applying rule 1 
when node 1 is considered both the initial and terminal node. Generating 
a second path requires backtracking to node 9, then continuing the sequence 
10-7-2-6-1. Note that the above graphical technique for listing all paths can 
be helpful when solving problems by hand. 

(3) 

b. A Detailed Description of the Path-Finding Algorithm 

Algorithm (Path-finding) ; This algorithm finds all paths between two 
nodes of a directed graph (without parallel edges) whose nodes are labelled 
1,2,...,N. The only modification necessary to adapt the algorithm to finding 

The format used to describe the path-finding algorithm follows the style of 

Knuth(^) , 
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all circuits thru node L is to set I < L where L and I are defined below. 

Notations: 


I: Initial path node 


L 

N 

E J 

R(J,M) 

C J 

P(V,W) 


u w 

F(K) 


Last path node 
Number of nodes in graph 
Number of branches leaving node J 
Routing table 

Column counter for the J row of the table R 
,th 


The V node in the node sequence of path W 
Number of nodes in path W 

A function used to test whether node K is repeated, and 
whether the last node is reached. 

PF1. (Preliminary) 

Set R(J,1), R(J,2) , . . . ,R(J,Ej) to the group of nodes of distance 

one from node J. When using the algorithm to find circuit, made 

the entries of each row decrease as M Increases. 

(--1 for M =» E t +1 and J=I 

Set R(J,M) < \ 0 f or M <3 e _+ l and J^I 


f 1 for K 83 ! 

Set F(K) < 1 0 for K=J and 

L-l for K“L 


Set Cj < 1 for J=l,2 N 

Set W < 1, V < 2, J < I, P(l,l) < I 

PF2 . (Find the next node) 

Set P(V,W) < R(J,Cj) 

PF3 . (Test R) 


IF R(J >Cj) 



< 0 stop; all paths have been found 

= 0 set F(J) < 0, go to step PF6 

> 0 go to step PF4 
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PF4. (Test F) 


r < 0 path completed; go to step PF7 
IF FtRCJ.C^H- 0 go to step PF5 

> 0 set Cj < Cj+1; 8° to step PF2 

PF5. (Prepare for next node) 


Set J < P(V,W) , F(J) < 1, V 


V+l, go to step PF2 


PF6 . (Back step) 

Set Cj <— 1, J <— P(V-2,W), Cj <— C +1, V <— V-l, go to step PF2 
PF7 . (Finish path) 

Set Cj <— Cj+1, P(K,W-1) <— P(K,W), K-1,2 l^-l, W <— W+l, go to 

step PF2. 
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II-4. Generating Nontouching Loops of Order Two or More 

Preliminary results from SNAP Indicate that of the following sub- 
programs, (1) finding a SFG, (2) coding and de-coding, (3) enumerating first 
order loops, and (4) finding all higher order nontouching loops, the last 
will generally require the most time unless the network contains many distinct 
symbols in which case subprogram (2) may dominate. It is therefore necessary 
to exercise considerable care in developing an algorithm for finding all 
orders of nontouching loops. 

In general, to find loop sets of all orders, some comparison between 
the node sequences of the different loops must be made. A brute force 
technique is simply to store all the node sequences of the first order loops 
and to find nontouching loops by direct comparison of the nodes contained 
in the loops. Of course, storage is also needed to indicate the loops 
contained in some of the higher order combinations, but this storage is 
necessary even in the more efficient techniques which follow. 

The above method is improved considerably if instead of directly comparing 
the nodes of loop A and loop B to determine if they touch, a function F(I) 
is defined as 


F(I) =* { 


1 

0 


Ie {nodes in loop A] 
otherwise 


and then tested as follows: 


If 



all Je {nodes in loop b} loops do not touch 
any Je {nodes in loop B] *=> loops touch 


For those computers which can accomodate the .AND. operation (or 
equivalent), the following coding technique reduces the number of logic 


Although the program was correct, the algorithm was incorrectly described in 
reference (3). 
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coraparisons needed to detennine if two loops touch to one and, perhaps 
what is even more important, requires only a single code be stored for 
each first order loop instead of the complete node sequence. As each first 
order loop is generated, it Is assigned an integer code whose binary repre- 
sentation shows the set of nodes in the loop. For example, if loop A 
contains the nodes {l 1 , 9,12, 8,5,6, l} and loop B contains the nodes {2,6), 
then the codes are 

A » <110110110001 ) 2 “ 3505 
B = (000000100010) 2 - 34 

To determine whether the two loops touch or not, the masking operation .AND. 
is used. Thus, 

(A) .AND. (B) - (000000100000) £ 0 

The result is not zero, indicating that loops A and B touch. 

Using the coding scheme the complete set of nontouching pairs of loops is 
found and stored in the one dimension array N. Let n » number of first 
order of loops. Then 

{N(l),N(2),...,NrP(l)],N[P(l)+l ],..., N[P(2)],N[P(2)+1] N[P(3)], 

N[P(3)+l),...,N[P(n -l)],N[p(n -l)+l], . . . ,N[P(n )]} 
is the complete set of nontouching pairs of loops, where the 

{N(l) ,N(2) , . . . ,N[P(1) ]) “ set of loops numbered higher than 1 which do not 

touch loop 1. 

{n[p( 1)+1 ],N[P(l)+2 ], . . . ,N[P(2) ]} “ set of loops numbered higher than 2 

l which do not touch loop 2. 

• 

{N[P(i-l)+l],N[P(i-l)+2t, ...,N[P(i)]} - set of loops nuebered higher than i 

l which do not touch loop i 

• 

{NTP(n -l)+l],N[P(n -l)+2 ], . . . ,N[P(n ) ]) = empty set because there are no 

loops numbered higher than n . 

Note that the array P is simply used to partition the array N such that the set 
{N[P(i-l)+l ],N[P(i-l)+2 ], . . . ,N[P(i) ]} does not touch loop i. 
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Gxaople: 

The 


(Consider the SFG of Fig. A) 
first order loops are 
loop node sequence 


1 

1-11-9-12-8-5 -6-1 

2 

1-11-9-10-7-2-6-1 

3 

2-6-2 

A 

2-3-2 

5 

A -5 -A 

6 

5-6-5 

7 

9-11-9 

8 

9-10-9 


To find the array of nontouching pairs P, SNAP codes 
proceeds to use the .AND. operator. The results are 


the above loops end 


N » {A, 5, 5, 7, 8, 5, 6, 7, 8, 7, 8, 7, 8} 

and P(l)-1, P(2)“2, P(3)°5, P(4)°9, P(5)»J1, P(6)-13, P(7)“13, P(8)=13. 

To systematically continue the process, an array S is created from which 
all higher order loop sets (2 or more) not touching loop 1 can be found. By 
Incrementing t from 1 to n, all higher order loops will then be enumerated. 
Let 

S(l,l)S(l,2)...s[l,U(l)] . . . S[1,J(1) 

s „ S(2,1)S(2,2) . . . . 8[2,U(2)+l]...S[2,J(2)],0,... 

S(K,1)S(K,2)...S[K,U(K)] . . . s[K,J(K)],0,... 


where S(l,l) - N[p(L-l)+l] 

S(1 , 2) - N[p(L-l)+2] 


s[i,J(i)] - N[p(L)] 
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and where the entires (loop lumbers) of row M(M > 2) are those loops in the set 

{sfM-l,U(M-l)+l ],S[M=l,U(M-l)+2] S[M-1, J(M-l) ]} 

which do not touch the loop S[M-1,U(M-1) ]. 

The arrows shorn in the array S given above are referred to os ’’pointers" . 

Note that U(J) indicates the position of "pointer" of the row. Example: 

ttl 

U(3) « 5 Beans the pointer of row 3 is currently located at the 5 column. 

The procedure for finding all higher order loop combinations is given in 
the following flow chart: 
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Generate row K+l of S as follows: Insert those 

loops of the set is[K,U(K)+l] S[K,J(K)]} 

that do not touch loop s[K,U(K)] into row K+l 
of S. Set J(K+1) ® number of these loops. When 
a new element, say S(K+1,X), is generated, the 
weight parameters corresponding to the symbol code, 
constant term, and power of s are stored (or when 
possible combined with other similar type terms) 
for the loop set loop 

loop i • loop S [l ,U( 1) ]•... -loop S(K+1,X) 


K“K+1 


Can we generate row K+l 
by incrementing the pointer 
of row K? i.e. Is U(K)<J(K)-1? 


Can we back up one row? 
i.e. Is K >1 
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Example: 

From the preceding example, 

N - (4, 5, 5, 7, 8, 5,6, 7, 8, 7, 8, 7, 8} 
and P(l) - 1, P(2) - 2, P(3) » 5, P(4) - 9, P(5) - 11 

P(6) « 13, P(7) «* 13, P(8) « 13 


Arrays N & P are more easily interpreted by setting up the following table: 

Table 3. 


loop J 

1 

2 

3 

4 

5 

6 

7 

8 


loops numbered higher than J that 
do not touch loop J 

4 

5 

5.7.8 

5. 6. 7. 8 

7.8 
7,8 


The sequence for producing the higher order loops is as follows ■ 
loop 1 

loops not touching loop 1 are inserted 
into first row of S (see Table 3) 

S array Output 


.if 

4 0... 
0 


no 3 r< * order loops 


loop 2 


S array 
r * 

5 0... 

0 


Output 

no 3 r< * order loops 
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f 


loop 3 

S array Output 

ri 

5 7 8-0 

0 


loop 5 does not touch loop 7 
or loop 8 (this is determined 
by comparing loop codes- -see 
V section II-4) 


loop 3 ‘loop 5 ’loop 7 
loop 3 ’loop 5 ‘loop 8 


loop 7 touches loop 8; thus, 
there is no 3 rd row. Further 
if the pointer of row 1 is 
incremented by 1, no new 2 nd 
row can be created. Thus, we 
^are done with loop 3 . 

loop 4 

S array Output 

J, 

5 6 7 8 

0 0 0 0 

* 


loop 5 does not touch loop 
v 7 or loop 8 

J 

5 6 7 8 

7 8 0 0 

0 0 0 0 
0 0 0 0 

loop 7 touches loop 8; thus, 
there is no 3 rd row. Incre- 
ment pointer of row 1 - 

v 


loop 4 ‘loop 5 ‘loop 7 
loop 4 ‘loop 5* loop 8 


5 7 8 0 

'l 

7 8 0 0 

0 0 0 0 

0 0 0 0 
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l 

5 6 7 8 

0 0 0 0 


loop 6 does not touch loop 
v7 or loop 8 

loop 4 ‘loop 6 ‘loop 7 
loop 4 ‘loop 6 ‘loop 8 


incrementing pointers give 
no additional third order 
v loops 

loop 5 

Output 

no 3 rc * order loops 


S array 
■ l 

17 8 0. 

0 


I 


loop 7 touches loop 8 


5 6 7 8 

7 8 0 0 
0 0 0 0 
0 0 0 0 


loop 6 

S array 

4 

7 8 0 0 

0 


^loop 7 touches loop 8 


Output 

no 3 rt * order loops 


/ 


loop 7 

S array 

[ 0 ] 

loop 8 

S array 

[0] 


Output 

no 3 r< * order loops 


Output 

no 3 r< * order loops 
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III. USER'S GUIDE 


I1I-1. Information Needed by User 


Program: SNAP (Symbolic Network Analysis Program) 


Purpose: To obtain the network functions 


V . V I I 

* out out out out 

V. ’I. ’ V. ’ ° r I 

in m in in 


as a 


ratio of two polynomials of the following type: 

(1) all network element values are represented by symbols (the symbols 


need not all 
V 

Examples: — 


be different) , 

2 

out _ s LRC 

in s 2 2LRC + s(L+R 2 C) + R 


V out _ ZYR 

V in 2ZYR + Z 4- R 2 Y + R 


( 2 ) 


(3) 


Description: Program SNAP is designed to handle lumped, linear, time invariant 

networks** containing the following type components: 

(1) two-terminal circuit elements -- resistance, inductance, and 
capacitance. 

(2) two-terminal networks described by an admittance or impedance 
parameter. 


some element values are specified numerically, some symbolically, 
V 


Example: 


s 2 R 


out _ 

V in s 2 2R + s(.5xl0 6 + 150R 2 ) + .75 x1Q 8 R 


all element values are given numerically, 


„ , out 

Example: — 2 


_ 4 8 

in 2s 4- 2x10 s + .75x10 


Refer to Appendix A at end of this section for a technique of handling 
multi-output functions. 

:ic 

See Appendix B for a brief list of additional limitations on the size and 
type of network allowed. 
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(3) all four types of controlled sources (Note: Mutual inductance, 

ideal transformers, gyrators, etc., can be modeled with elements 
in (1) and (3)) 

(A) one independent source; see Appendix A for a technique of handling 
multi-input networks. 

Network Data Required: After the network components have been modeled by the 

type elements allowed, the branches and nodes are to be numbered consecutively 
starting with 1 and reference directions for each branch current are to be 
chosen. The following gives the sequence of data cards needed to describe 
the network. 

CARD 1 

Columns Contents 

1-72 Title card (all 72 columns are reproduced in output) 

CARD 2 

Columns Contents 


1-5 

/ right . 
^adjusted' 

6-10 

right 

'‘adjusted'' 


Number of nodes in the network 


Number of branches in the network 


The following three entries are optional. 

11-15 

, right . Number base of symbol codes (automatically set to 8 if left 

''adjusted'^ blank) 


21 


22 


1 if a description of the SFG is to be listed, blank otherwise 

1 if all loops (circuits) in the SFG are to be listed (node 
sequence), blank otherwise 



CARD 3 


Columns 


Contents 


1-5 

/ right . 
'‘ad justed ; 

6-10 
/ right 
adjusted' 


Network branch number of source 


Network branch number associated with output (leave blank 
if output is a voltage across more than one branch) 


11-15 
right j 
'“adjusted'^ 


Node number corresponding to the positive output voltage 
terminal (these columns can be left blank if columns 6-10 
are not blank) 


16-20 
, right 
'’adjusted' 


Node number corresponding to the negative output voltage 
terminal (these columns can be left blank if columns 6-10 
are not blank) 


CARDS 4 thru (b+3) 

(b = number of network branches) 

Note 1: Each card describes one network branch (element). 

Note 2: If output is a voltage (current) associated with a particular branch, 

then the data card describing this branch should be entered first 
(last) among the branch data cards (cards 4 thru (b+3)) to insure 


that this branch will be chosen as part of the tree (cotree). 
Note 3: When a large number of branches share one common terminal, it is 


better to place these branches first starting with card 4 (card 5 
if note 2 applies). The reason is given in Appendix C at the end 


of this section. 


Columns 


1-2 


( 


left 

adjusted 


) 


Element type.; E 
1 
G 
R 
L 
C 
Z 
Y 
CC 
CV 
VC 
VV 


Contents 

voltage source 
current source 
conductance 
resistance 
indue tance 
capacitance 
impedance 
admittance 

current controlled current source 
current controlled voltage source 
voltage controlled current source 
voltage controlled voltage source 
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Continued 

Columns 

3-5 

right 

''adjusted'* 


Contents 

Element number--all elements of the network must be assigned 
a distinct number (positive integer). For greatest efficiency, 
the numbering should be consecutive. 


6-10 

right 

^adjusted' 


Initial node--this is relative to the arbitrarily chosen 
current direction. 


11-15 
, right 
'adjusted' 


Terminal node--this is relative to the arbitrarily chosen 
current direction. 


17-19 
/ right . 
''adjusted' 


Element symbol--the element's value, if not specified, is 
represented by this symbol. 


20 


21-32 

right 

'“adjusted' 


Equal sign (=) if element is to be assigned a value. Leave 
blank if element value is to be represented in symbolic form. 

Element value (if known) --Format is E12.5. Units should be 
compatible with element type as specified in columns 1-2; for 
example, R is expressed in ohms, G in mhos. 


33-35 

( right . 
''adjusted'* 


If element is a dependent source, enter the element number of 
its control. 
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DATA DECK 


R 

11 

5 

1 

R3=> 

10. E3 


R 

10 

8 

5 

RE=25. 



R 

9 

8 

5 

RE=85. 



CC 

8 

4 

8 

A2 


9 

CC 

7 

3 

8 

Al 


10 

R 

6 

3 

4 

R2« 

1S.E3 


C 

4 

1 

4 

C 



R 

5 

1 

4 

Rl = 

5.E3 


C 

3 

1 

3 

C 



R 

2 

1 

3 

Rl= 

5. E3 


E 

1 

8 

I 





1 

6 






5 

It 


11 




\ FFERENT 1 AL flMPL ! F 


n n ft n n n ft n n n 


(i 1 1 n ft ii n 


D 0 8 0 0 5 1- 8 8 0 e (1 8 C B 0 5 0 I 0 0 S C 3 D 5 0 !! 9 e ? C C 0 0 B D P 6 0 0 C P e C e 0 (! 0 0 0 8 0 C 0 0 P 0 0 0 0 G 0 0 0 9 C C D 0 D C 6 0 ; ? 0 

I j 3 « 5 t 7 e J hU! I’ (3 ><!■■ ;•>’>? 1! p:' 7:-’t S7 7*'K 7>-‘C 51 St. ’5 }! 75 2S :? it ’J C 4U/ <3 </ 15 :5« 75 « 5: 5' aSHSSiK 57 ’C'SHUI f.n 13 11 7? h,lu 'll ■■ j 

1 1 ] 1 ] 1 H 1 1 1 1 11 1 1 i i 11 ? i 1 ? 1 1 1 1 1 1 1 I 1 11 1 1 1 11 I i i 1 1 1 1 1 I 1 ! I 1 I l 1 1 1 I !1 I H I 1 1 11 1 1 1 11 M ! I i 


TABLE 1. Input Data as Reproduced in Program Output. 


<mmh»o®<h>««dIfFERF.NTI AL AMPlIF I eR*#****#*** 


NUMBER OF NODESs 5 

NUMBER OF BRANCHES 5 11 

ELEMENT NUMBER OF SOURCE 5 1 

ELEMENT NUMBER ASSOCIATED WITH OUTPUT * b 

BASE FOR SYMBOL CODES 5 p 


NETWORK 


ELEMENT 

ELEMENT 

INITIAL 

terminal 

ELEMENT element elemfnt Nn 

type 

NUMBER 

NODE 

node 

symbol VALUE of 

CONTROL 

E 

1 

2 

l 

-o. 

-0 

R 

2 

1 

3 

Rl 5 5.00000E*03 

-0 

C 

3 

1 

3 

C -0. 

-0 

R 

5 

1 

4 

R 1 5 5.00000F + 03 

-0 

C 

4 

1 

4 

c -o. 

-0 

R 

6 

3 

4 

R2= 1.80000F*04 

" 0 

CC 

7 

3 

2 

Al -0. 

1 0 

CC 

8 

4 

2 

A 2 “ 0 « 

9 

R 

9 

2 

5 

RE 5 ?.50000F*01 

-0 

R 

10 

2 

5 

RE 5 2.s0000E*01 

•» 0 

R 

11 

5 

1 

R3s 1,OOOOOE<-04 

a 0 




TREE 

selected 


E 

1 

2 

1 

-o* 

“0 

R 

2 

1 

3 

R 1 5 5. 00000E*Q3 

“0 

R 

5 

1 

4 

R 1 5 S, 00000F*03 

-0 

R 

9 

2 

5 

RE 5 P.SOOOOE+Ol 

" 0 
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TABLE 2. Program Output Information Showing Signal-Flow Graph, 
Circuits, and Execution Times. 


SFG 


NITIAL 

TERMINAL 

EXPONENT 

BRANCH 1 

BRANCH 

1 IF SYMBOL 

1 IF SYMBOL 

NODE 

node 

OF S 

Value : 

symbol 

IS inverted 

IS USEO 

6 

1 

0 

-1 o00000E>00 

FB 

0 

1 

2 

3 

1 

1 , 00000£*00 

C 

0 

1 

3 

2 

0 

~5«OOOOOE*03 

R 1 

0 

0 

5 

4 

1 

1 ,oooooe*oo 

C 

0 

1 

A 

5 

0 

“5 o 00000E+03 

R 1 

0 

0 

2 

6 

0 

-6,66667E“05 

R2 

1 

0 

6 

2 

0 

5* OOOOOE+03 

Rl 

0 

0 

5 

6 

0 

6.66667E-05 

R2 

1 

>0 

6 

5 

0 

-5« 00000E + 03 

Rl 

0 

0 

10 

7 

0 

1.00000E*00 

41 

0 

1 

7 

2 

0 

5.00000e*03 

Rl 

0 

0 

9 

12 

0 

4.00000E-02 

RE 

1 

0 

12 

8 

0 

1 • OOOOOE* 00 

A? 

0 

1 

8 

5 

0 

5. 00000E*03 

Rl 

u 

0 

9 

10 

0 

4 « 0 0 0 0 0E“ 0 2 

PE 

1 

0 

10 

9 

0 

-2.50000E*01 

RE 

0 

0 

9 

11 

0 

-1, 0O000E-04 

R 3 

1 

0 

11 

9 

0 

2.5000 OE *01 

RE 

0 

0 

1 

11 

0 

1.00000E-04 

R3 

1 

0 


TIME for FORMULATING SIGNAL 

FLO* GRAPH IN SECONDS *25? 

CIRCUITS 


NO* 

NODE 

list 



1 

1 

11 

9 

12 

8 5 6 1 

2 

1 

11 

9 

10 

7 2 6 1 

3 

2 

6 

2 



4 

2 

3 

2 



5 

4 

5 

4 



6 

5 

6 

5 



7 

9 

11 

9 



8 

9 

10 

9 



time for 

FINDING 




8 FIRST ORDFR LOOPS 

in seconds 


0 

046 


time for 

finding 




1 9 SETS OF 


NONTOUCHING LOOPS? IN SECONDS .027 


TIME FOR DECODING SYMBOLS IN SECONDS 


.124 
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TABLE 3. Network Transfer Function and Total Execution Time. 




COLUMN 

1 

2 

3 

4 


NUMERATOR POLYNOMIAL 

- (j, 33333 xio'^-t- .!&&&(, 7s C)(A3~A l) 

SYMBOL FOR GIVEN COLUMN 
A2 / l 

Al / 1 

C A2 /I 

C Al / 1 


POWER 
OF S 

COLUMN 1 

0 3 ® 33333E** 05 

1 0, 

2 0 . 


CONSTANT COEFSe 
COLUMN 2 
“3 « 33333E“05 

0, 

0. 


in the polynomial 


COLUMN 3 

0. 

1 .66667E-01 

0. 


COLUMN 4 

0. 

- 1 « 66667E- 0 1 

0. 




denominator POLYNOMIAL 

— 3.33 7 jr -h £ 7X 10*~ 5 (L f S. OO &3SK! O^S^Cf" 


COLUMN 

SYMBOL for GIVEN 

COLJMN 

1 

1 / 

1 

2 

c / 

1 

3 

C*<»2 / 

1 


POWER 
OF S 

COLUMN 1 

CONSTANT COEFS, IN THE POLYNOMIAL 
COLUMN 2 COLUMN 3 C 

0 

3»33750E*00 

0 • 0 * 

1 

0. 

2.F.7000E + 04 0, 

2 

0. 

0. 5, 00625E ♦ 07 


EXECUTION TIME IN SECONDS, ,509 


AUGUST 1970 VERSION OF SNAP 
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APPENDIX A 

A Sorting Technique for Handling Multi-Input, Multi-Output Networks. 
Multi-Inputs 

Program SNAP (the August 1970 revision) permits only one independent source 

branch. However, networks containing more than one source can easily be 

handled with the following technique. Let i = l,2,...,n represent a set of 

n independent sources, either voltage or current. Assign as the permitted 

independent source and make W„ , W„,...,W dependent sources which are dependent 
' IS n — 

on Wj with proportionality factors 

W 2 W 3 W n 

k 2 ' ■ k 3 * k „ “ respectively 


Only the numerator polynomial in the output will contain these parameters thus 
permitting the user to easily put the output function into the form 


!out 

W 1 

where A and P^ 
be written 


Pi t r 2 k 2 -e ^ -k . . + P n k n 

A 

i = l,2,...,n are polynomials. 


The output function can then 


W 


P..W. + P„W„ +. . .+ P w 
11 2 2 n n 


( 1 ) 


out A 

Although at present SNAP does not give the output function in the form of 

Eq. (1) directly, only a few program modifications are necessary to effect such 

a result. For example, the program could internally create a new input node, 

I , of the SFG and then make each independent source, W. , dependent on I 
new r l r new 


with weight P^ as shown in Fig. 1 below. 
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Figure 1 


Multi-Outputs 


The following technique can be used to obtain more than one output function 
in a single computer run: Augment the original network by appending one end 
of a series connection of dependent voltage sources to the given network such 


(a) to each branch current, Ij, desired as an output, there corresponds 

a dependent voltage source which depends on 1 ^ and has symbolic weight 


and (b) to each voltage desired as an output, there corresponds a set of the 
dependent voltage sources each dependent upon a voltage across one 
of the branches in the path between nodes A and B and all having 
symbolic weight V^g. 

By specifying the output to be the voltage across the entire series connection of 
dependent voltage sources, outputs 1^ and V^g will be those output terms which 
contain 1^ and V^g respectively. Only a few modifications of the present 
version of SNAP would be necessary to have the program internally perform the 
network augmentation described above (at present, the user must do the augmenting). 

As an example, Fig. 2 illustrates the network augmentation needed to find 
the voltage V ^ and current 1^ for the given bridge network in one computer run. 




A Brief List of Limitations on the Size and Type of Network Allowed 
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APPENDIX C 


Selecting a "Good" Tree 

The network tree used to generate the SFG has a very significant effect 
on the number of loops and higher order loops present in the SFG. The 
loop enumeration and evaluation, in turn, often determines the time and 
storage needed by a computer to solve a given network. The ladder network 
of Figure 1 together with Table 1 illustrated the interrelationship between 
the tree selected, the number of loops (all orders), computer execution 
time, and computer storage. 
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Table 1 

Tree Branches Number of Number of Time required to find V /I 

loops higher order ° b 

loops 


Star tree: 1,3, 5, 7, 9, 

11,13,15,17 

17 

2567 

1.55 

seconds 

1,2,4,7,9,11,13,15,17 

38 

8096 

3.93 

seconds 

1,2,4,6,9,11,13,15,17 

117 

19719 

9.42 

seconds 

1,2,4,6,8,11,13,15,17 

476 

(storage for first 
order exceeded) 


-- 


Unfortunately, choosing the "best" tree, that is, a tree which will 
minimize the number of loop combinations of all orders is a very involved 
process. See reference 2 and 7 for a detailed discussion of this problem. 
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< 

For most networks, however, a tree that will result in a reasonable amount 

of execution time and computer storage can be selected by applying one 

of the following rules (rule 2 results in a better tree than rule 1) 

Rule 1 : Select a tree in which as many branches as possible form a star, 

that is, the branches share a common node. Modify this tree, if necessary, 

to include any branch which has two or more branches in parallel with it. 

Rule 2 : Let be some tree (not necessarily the best) of the network 

graph. For each link X. of the graph, define as the number of tree 

i 

branches which form a circuit with l . . Then form the sum 

l 

L 

S^ = £ Bp where L = number of links in the graph having 

k i=l i 

as a tree. 

Select that tree, say T. , which satisfies the inequality 

St < S k = 1,2,...,N where N = number of trees 
j T k 

The example given in section 1II-1 uses a tree, T., having S^ =11. 

J J 

The tree generated internally by program SNAP includes all voltage 
sources together with those passive branches read in first (starting with 
input data card 4) which complete the tree. Thus, to have SNAP select 
the tree that has been chosen by the user, it is necessary that the user's 
tree include all voltage sources and that all its passive branches be 
listed first starting with input data card 4. 
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III-2. Modifying the Dimension of Arrays 

In order to make SNAP applicable to many different type networks, a flexible 
yet simple procedure is needed for modifying the dimension of the arrays. For 
example, storage requirements for networks containing many symbols will be 
determined by the number of symbols, symbol codes, etc., whereas the storage 
needed for networks having no symbols will be determined by the number of loops, 
nontouching loops of all orders, and related network characteristics. Because 
it is not possible to determine apriori reasonable bounds for all the network 
characteristics, error diagnostics have been built into the program to inform 
the user as to which arrays have been inadequately dimensioned. As a result, 
the technique for adjusting the array dimension, in SNAP can be outlined as 
follows : 

(1) Check that those network characteristics which can be determined before 
running the program are within the specified limits. These limits are 
listed following the dimension statements of the main program for 
convenient reference. 

(2) Run the program. If an array dimension is exceeded an error message will 
result which specifies the network characteristic involved. For example, 
if the SFG of a given network has an excessive number of circuits, the 
message "No. of circuits exceeds limit--increase dimensions containing 
NPAC" will result. The definition of NPAC (number of paths and circuits) 
are found immediately following the array dimensions in the main program. 

It is important to point out that a computer run may continue to completion 
even if the dimension of some arrays have been exceeded (an error message is 
still given, however). In this situation, the results cannot be considered 


reliable , 
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(3) Once it has been ascertained by (1) and (2) that dimension modifications 

are in order, refer to the next few pages to determine the arrays associated 
with the network characteristics of interest. Increase the dimension of 
all the arrays indicated by say 207„ (several runs may be necessary to 
achieve adequate program dimensions). Then update the value of the 
parameter (NPAC, for example) corresponding to the network characteristic 
involved. This parameter is used throughout the program (as limits on 
DO loops etc.) thereby making it unnecessary to do any additional program 
modifications. 


NBN ■= Number of Network Branches (Presently 35] 




PROGRAM 

MAIN 

IG(NBN) , 
SMBOL(NBN) , 
I FLOW (NBN) , 
LT(NBN) , 

KODES (NBN) , 
KONC(NBN) , 
N(NBN,NBN) , 
NP(NBN) , 

KODE(NBN,NBN) 
IXPO(NBN,NBN) 
CONS (NBN, NBN) 




SUBROUTINE SFG 

J ROW (NBN) , 
NP(NBN) , 
IVV(NBN) , 
NUML(NBN) , 
ICV(NBN) , 
INTREE (NBN) , 
LINC(NBN) , 

TYPB(NBN) , 
JB(NBN) , 

LB (NBN) , 
MSYM(NBN) , 
IQUAL(NBN) , 
VAL(NBN) , 
SYM(NBN) , 

IQUALX(NBN) , 
VALX(NBN) , 
NUMLX(NBN) , 
INTRE (NBN) , 
NOTREE (NBN) , 
TYPX(NBN) , 
NUMX(NBN) , 

JBX(NBN) 
LBX(NBN) 
IB(NBN , NBN) 
NS (NBN, NBN) 
NF(NBN,NBN) 



SUBROUTINE FTREE 

TYPX(NBN) , 
JBX(NBN) , 
LBX(NBN) , 

INTRE (NBN) , 
NO TREE (NBN) , 
NP(NBN) 

NF (NBN, NBN) 



SUBROUTINE TREP 


JX(NBN) , 
NP(NBN) , 


JMEM(NBN) , NF (NBN, NBN) 
KMEM(NBN) 
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NBG = Number 


NFIRST(NBG) , 
NLAST(NBG) , 
IXPON(NBG) , 
WEIGHT(NBG) , 


NFIRST(NBG) , 
NLAST(NBG) , 
IXPON(NBG) , 
WEIGHT(NBG) , 


NPAC = Number 


CONST (NPAC) , 
KODET(NPAC), 
IXPOT(NPAC) , 


f Branches in SFG (Presently 100) 

PROGRAM MAIN 


SYMBUL(NBG) , NEST(NBG) 
MIX(NBG) , TYPE(NBG) 
CVAL(NBG) , 

KONSO(NBG) , 


SUBROUTINE SFG 


MAPY(NBG), SYMBUL(NBG) 
KONSO(NBG) , MIX(NBG) 
NEST(NBG) , CVAL(NBG) 
TYPE (NBG) , 


of Paths Plus Circuits (Presently 300) 

PROGRAM MAIN 

MAPO(NPAC) , JAC(NPAC) 

NOCTOT(NPAC), NPCODE(NPAC) 

NUP(NPAC) , 
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NTO = Number 

of Terms in Output 

(Presently 150) 

NA(NTO) , 

POLYU(NEXPS.NTO) , 

PROGRAM MAIN 
SEMPON(NTO,NSPT/2) 

NB(NTO) , 

POLY (NEXPS, NTO) , 

SEMPOD (NTO , NSPT/ 2 

KSORT(NTO) , 

SIMB0N(NT0,NSPT/2) , 

ITOP(NTO) , 

SIMBOD(NTO,NSPT/2) 


KSORT(NTO) , 

POLY (NEXPS , NTO) 

SUBROUTINE ARRAY 

I TOP (NTO) 


SUBROUTINE DECODE 


NSPT = Number of Symbols per term in Output (Presently 20) 

PROGRAM MAIN 

KONS(NSPT) , KODF(NSPT) , SEMP0N(NT0 ,NSPT/ 2) 

KODI(NSPT), SIMBON(NTO,NSPT/2), SEMP0D(NT0,NSPT/2) 

SEMBOL(NSPT) , SIMB0LD(NT0,NSPT/2) , 

SUBROUTINE FTREE 


KCOL(NSPT) 


SUBROUTINE DECODE 


SEMBOL(NSPT) , KODF(NSPT) , 


KODI(NSPT) 
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NEXPS = Number of Different Powers of s (Presently 15) 

PROGRAM MAIN 

MSORT(NEXPS) , POLYU(NEXPS ,NTO) ,POLY(NEXPS,NTO) 

SUBROUTINE ARRAY 

MSORT(NEXPS) ,POLY(NEXPS,NTO) 

NRI = Maximum Number of Nontouching Loops (Presently 15) 

PROGRAM MAIN 


ISET(NRI , NCI) 

NCI = Maximum Number of Loops Not Touching any Given Loop (Presently 100) 

PROGRAM MAIN 


ISET(NRI,NCI) 

NEON = Number of Nontouching Paris of Loops (Presently 1200) 

PROGRAM MAIN 


NOTCH(NEON) 

NRS = Number of Repeated Symbols (Presently 9) 

PROGRAM MAIN 


STAR(NRS) 
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IV. PROGRAMMER'S GUIDE 
IV- 1. Definitions 


CONS ( J»l> sWEIGT ( I ) FOR BRANCH I OF THE SFG WHERE 
J= FIRST ( i ) * L=NLAS I (I ) 

CONST (I)=C0mPOSt1E CONSTANT ASSOCIATED with circuit t. it TF found hy 
TAKING THE PRODUCT OF THE CONSTANT VALUES OF EVERY SFG BRANCH 
IN CIRCUIT I 

CVAL(NUMC) sVALX (LINK) WHERE NUMC*NUMX (LINK) 

(USED ONLy FOR NETWORK BRANCHES NOT IN THE TREF) 
IB(LF»JF)=I8 (JF,LF)sNUHC WHEhE JFrJB(NUMC) AND IML«(NUMC) 

AMO NUNC IS A NF T WORK I REF BRANCH NUMBER (ASSIGNED P,Y USFR) 
IFI.Ow(K)=A FLAG, FOR THE PURPOSE of CHECKING WHFThFR NOOF K TS PFPFATFO 
AND WHETHER THE LAS I NODE IS REACHED 
TG(L)=SYMBOL COot ASSIGNED To THE SFG BRANCHES RAVING 
TFRMJNAL NnUF L 

INTRF(K)aI, THL I-Th NfTWORK BRANCH TN THE DATA BRANCH LIST TS 
CHOSEN AS The K-Th BRANCH OF The network tree 
TNTRFF (MlJMC) =1 tF THE NETwORk BRANCH NUMBERED NUmC BY THE USER IS 
SELECTED FUR THE TREE » 0 OTHERWISE 
IUUAL (NUMC) sIuUaLX ( I ) WHERE NUMC=NUMX(I) 

(USED ONLY FOR NETWORK TREE BRANCHFS) 

IOUAL.X ( I ) sFQUAl STGN(s) TF UTH NETWORK BRANCH IN THF OaT A BRANCH 
LIST HAS a NUMERICAL VaLUE. LEFT HLANK IF I-Th Branch ts 
TO HE REPRESENTED BY A SYMBOL 

iset ( j» i ) sThf integer array which together with the array notch r»*' 

BE USED Tn FIND AlL SE 1 S OF NONTOUCHING LOuPs OF ORUFR GPFATFR 
THAN 2 

I Top ( jc) s i if ThE TfcRMs in column jc of the array poi.y belong to 
the numerator of the output transfer function* o tf thfy 

BELONG 10 The DENomINAIOR 

TW (M) SNFTW0RK RpANCH NUmRER OF ThE M-Th VOL'fAGF CONTWOLEO 

voltage Source in i he data branch list 

I XPO ( J*L) =1 XPON ( I ) FOR BRANCH I OF T HE SFG WHERE 
JsNFIHST ( i ) * LsNLAST ( I ) 

TXP0N(I)=FXP0NEm1 OF S ASSOCIATED WITH THE VALUE OF THE SFG HPANCH I 
TXP0T(I)sC0mPOSt If EXPONENT OF s FOR CIRCUIT I. IT TS FOUND BY ADDING 
THE 5 POWERS ASSOCIATED WITH Each branch In CIRCUIT T 
UAC ( J) snUmbeh Of NONZERO entries IN row j of isft 

JH (NUMC) sjHX ( I > WHERE NlJMC = Ni jM X ( I ) 

(USE 0 ONLY FOR NETWORK TR r F. BRANCHES) 
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JflX ( I ) =INITIAL MODE OF THE I-TH NFTWORK BRANCH IN THF.OATA 
BRANCH L I =>1 

JMFM(I)sTHE ROW of THE ROUTING MATRIX FROM WhlCH the i-th nodf 
IN THE path SEQUENCE was taken 

JR0W(LF)=THE NU.ibER OF NON-ZtHO ENTRIES IN ROW |_F OF THE ARRAY MF 
JX(l4l)sNP(I) 

kbasissnumher b a se of the symbol codfs. THAr is* thf sfg contains koo 
DISTINCT symbols* SEM80L (K) «K=1 *KOO* not INCLUDING THF ear 
variable s* WHERE SEMdoL(K) IS ASSIGNED THE CODE kBASTS*«K 
KHOLsCOUNTER UStD TO FIND THF NUMBER OF LOOPS OF ORDFR ? OR 1 GRFATFR 
KlKeA ROW COUNTER Oh THE MATRIX POLY 

KMEM(I)=THE COLUMN of Thf ROUTING MATRIX from which thf i-th nodf 
IN THE PA t H SEQUENCE WAS TAKEN 

KOOE ( J*L) -CODE REPRESENTING | HE SYMBOL OF IHE SFG BRANCH HAVING 
J AS AN INITIAL NODE AND L AS THE TERMINAL NODF 
KODES ( J) =2"* ( J-l ) WHEREJ IS a nODF OF THF SFG 

KODET ( I ) ^COMPOSITE CODE ASSOCIATED WITH CIRCUIT I. THIS CODE PEPRFSF N 
THF SET OF SYMBOLS LORrESdOMDInG TO I HE SET OF SFG BRANCHES 
CONTAINED IN CIRCUIT I 

KOI)F(NZ) IS THE MULTIPLICITY of the SYMBOL CORRESPONDING To THE 
CUOE KOOT(NZ) 

KOO I (NZ) *NZ*1 *2» . . . « IZ IS THE SET OF INDIVIDUAL SYMBOL CODES THAT 
MAKE UP ThE COMPOSITE CODE KSORT(JZ) 

KONC ( J) =COLUMN cCUNlER FOR ROW J OF THE ROUTING MATRTX N(J,K) 

konsckozy)si if the symbol having code kozy is not am inverse symbol a 

0 IF THE SET OF SYMBOLS CORHE5PONOT n 6 TO ThE COMPQSITF CODE 
KSORT(J) belongs to The DENOMINATOR polynomial 
konso { I ) s i if symbol of the sfs branch i=symbul(i), n if symbol of 
the sfg Branch i=i/symhul(I> 

K SORT ( K ) sTHE CODE ASSIGNED To COLUMN K OF THE MATRIX POLY 
LB (NUMC) sLBX ( 1 ) WHERE NUmC=NUMX(I) 

(USED ONLY EOK NETWORK TRrE BRANCHES) 

LBX ( I ) sTFRMINAL NODE OF THE t-TH NETWORK BRANCH IN THE DATA 
BRANCH L I S I 

LlLsA COLUMN COUNTER OF THE MATRIX POLY 

LINC(NUMC)31 if t he network branch numbered numc By the USFP is mot 
IN THE IRFfc* o OTHERWISE 
LIST=NUMHFr OF DIRECTED BRANCHES in the SFG 

LlSTC=l IF ALL rlRCUlTS OF T hE sFG ARF TO BE LISTED IN THF PRINTOUT, 

0 otherwise 

LlSTGsl TF SFG 1 NF ORMA T I ON ( dW ANJCH S YmROL S * wE i CRTS ETC.) ARf TO Rf 
LISTED IN THE PRINTOUT, 0 OTHERWISE 
LlSTPsl IF ALL PATHS FROM NODE NIN TO NODF. NUUT aRE TO RE I ISTFD IN 
THE PRTNTpUT, 0 OTHERWISE 

LT ( J) cNUmRER OF POSITIVE ENTRIES IN ROW J OF N(J,K) 
MAPO(NlP)sNOCT0T(NlP)»NncT0T{NlP-l) WHICH EQUALS THE NUMBER OF 
LOOPS NUT TOUCHING LOOP NIP 

MlX(I)=MAPPING OF THE SFG BRANCH LIST INTO A LIST SATISFYING ONF OF 

the following conditions 

NFIRST ( J) .GT.NFIRST (K) FOR J*GT,K 

OR NFIRST ( J) =NFIRST (K) , Nl AS H J) . LT , NLAST ( K) FOR J.GT.K 
M$ORT(K)=TH£ EXPONENT OF S ASSIGNED TO ROW K OF 1 Hp MATRIX POI Y 
N(J»K)» WHERE K-l ,2* . . • ,LT ( J) , IS THE TERMINAL NUDE OF . SFG BRANCH 
HAVING J aS ITS INITIAL NODE* THF VALUE Of- EACH NONZERO EMT»Y 
IN A GIVEm ROW IS MADE TO DECREASE AS K INCREASES. THF ADDITION 
ENTRY N (N t N,LT (NIN) ♦ 1 ) s-l IS ALSO MADE 
nA(J)sNUMBER OF SYMBOLS (NOT COUNTING INVERSE SYMBOLS) IN THF 
CODE KSORtU) 

NB ( J) sNUMBER UF INVERSE SYMBOLS IN THF COOL KSORl(J) 

NCIR=1 IF CIRCUllS ARE TO HE FOUND, AND 0 IF CIRCUITS ARE MOT TO 
BE FOUND 
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nEsT(T)= 1 IF THp sfg branch I CONTAINS A SYMBOL in addition to thf 
laplace variable s* o if the sfg rwanch i contains no symbol 

EXCEPT POSSIBLY FOR T Ht LAPLACE VARIABLE S 
NF (LF* JBOJ) =ROU T Ing TABLE FOk THE NETWORK COMPOSED ONLY OF 
BRANCHES wtLONGING TO I HE TREE 

nf i R" l if paths are to be found (nf ir set to i if listp*d* ano o tf 
paths are not to be found 

NFIRST (I)=lNiri *L NODE OF T Ht DIRECTED SFG (BRANCH I 

NlNsNETWORK branch NUMBER of THF SOURCE. this RtCOMFS THF SOURCE w OOF 
OF THE SFg 

NLAST ( T )= terminal NUDE of IHt DIRECTED ST G BRANCH T 
nohsnumbkr of branches tn neiwqrk 
NOOsNUMBF.H OF NoOF.S in NETWORK 

NODAsO unless OUTPUT is a VOLTAGE taken ACROSS MORE than onf. NETWORK 
ELEMENT. IN THIS CASE IT DEIGNATES THE POSITIVE TERMINAL OF 
THE OUTPUT VOLTAGE 

noobsO unless output is a voltage taken across more than one nftwork 
ELEMENT. IN this CASK It DESIGNATES THE NEGATIVE TFRmInAI 
OF THE OUTPUT VOLTAGE 
NOLsNlJMBER OE CIRCUITS (LOOPS) 

NOPsNUMBER OF PaThS F ROM NODE NlN TO NODF NOUT In THF SFG 
nout=netwohk branch number associated with the output (voltagf across 

OF CURRENT THRE) . THIS BECOMES THF SFG NODE CORRE SPONO T NG TO I 
OUTPUT VARIABLE 

NOTCH ( NOC ) ANO mOCTOT (K) . CONSIDER THE INTEGER $ET 

(I)s(l*2*.«**N2) WHERE N2SNUM8ER OF NONTOUCHTNG PaI»S OF LOOPS. 
NOW CONS I uE R THE FOLLOWING SUBSET OF (1). 

S ( I ) = ( NOC f O I (K-l) ♦1*NOCTOT (K-l J *2, . . . « (mOCTOT (K) ) WHEOF 
NOCTOT <rj) =0. IHEn I he SET (NOTCH(J>* J IN S(D) JS THF SET OF 
LOOPS That do not touch loop k 

NOTREE ( T ) — 1 IE r HE I°TH NfTWoRK BRANCH IN THE DATA LIST I s CHOSEN 

for the rptEi o otherwise 

nP ( I ) hTHF NODE SEQUENCE of a path between node NlN AMO NODE NDDT 

OF THE SF G » IF NINsNOuT THIS IS THE NODE SEQUENCE FDR A CIRCUIT 
OPCODE (K) sCOMPOSl TE CODE USED TO IDENTIFY CIRCUIT I. FOUND HY SUMMING 
THE CODES, KOUtS ( J) * AlLOTED Tq EACH NODE * J, TN THF CIRCUIT 
NIJML (NUMC) sNUMLA ( I ) WHERE NUmCskjUmX ( I ) 

(USED ONLY FOR NETWORK TREE BRANCHFS) 

NS ( LE » JF ) a 1 IF THE. NETWORK TREE BRANCH T B ( L F ♦ J F ) HAS INTIIAL NODF 
LF AND TERMINAL node jf and EQUALS -1 IF The nftwork trfe 
BRANCH HAs INITIAL NOOt JF AND TERMINAL NOuE LF 
N'JMLX(I). IF I»1 h NETWORK BRANCH IN THE DATA BRANCH LIST TS A 

DEPENDENT SOURCE* THIS ARRAY EQUALS THE NETWORK BRANCH NUMBER 

assigned to its control 

NUMX(I)=THE nE T t.iDRK BRANCH NUMBER ASSIGNED BY THt USED TO THF T-TH 
NETWORK HpAwCH In THE U A T A BRANCH LIST 
NUp(J) DESIGNATES THE LOOP* ISFT ( J*NUP( J) ) * OF ROW J WHICH IS MOT 
TOUCHED by THE LOOPS EnTEREO In ROW J*i OF IsET. 

ROLY (K »L) sMATRI X OF CONSTANTS WHERE EACH ENTRY TS ASSOCTaTFD WITH A 
TERM in T hE NUMERATOR OR DENOMINATOR OUTPUT POLYNOMlAI 
HAVING THr S POWER OF K AND THF SYMBOL CODE ASSIGNED TO COLUMN L 
POLYIJ (K,L) =MATNiX OF CONSTANlS WHERE EACH ENTRY IS ASSOCIATED W T TH A 

TERM IN THt NUMFRATOR OF tHF OUTPUT POLYNOMIAL HAVING THE S PDWF 
OE K AND THL SYMBOL COl.)E ASSIGNED TO COLUMN L 

skmhol ( ko) bsymBaL corresponding To The coot kbasis** iko-i > 

Sb MPON ( J1 » U2 ) »J?=1 *2*...»NA(J1) * AND SEmPOIJ ( J 1 * J 3 ) , J3 = l *2* . . . *NB Ml ) 
ARE RESHEc T TVtLY THE MULTIPLICITY OF THF. SYMBOl S 
SIMHON ( Jl ,J2) »J2 s 1*.«.nA(J]) t AND STMHOD < Jl * J3) » J3=l * . . . * NH ( .1] ) 
SIMBON ( Jl » J2) * J2= l 9 2* • O . * NA (Jl ) , aND S I MHOD ( J 1 , J3 ) , J3= 1 , 2 * . . . , NR ( J 1 ) 
ARE RESPECTIVELY THE SYMBOLS AND INVERSE SYMhOI S CORRESPOND 1 NG 
TO THE SYMBOL CODE KSOpT(jl) 
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SMBOL <K) sSYMBUL ( 1 > FOR THE SFG BRANCH I WHERE IsmIX(K) 

STAR { I ) »** I This array is generated from data statements and is 

US D IN FORMING THE ARRAYS SEMPON AND SEMPOD 
SYM(NUMC) sSYMX ( T> WHERE NUMCsNUMX ( I > 

(USED ONLY FOR NETWORK TREE BRANCHES) 

SYMRUL(I)=SYMrtOL ASSOCIATED WITH THf VALUE OF THE SFG BRANCH T 
SYMXtI)sSYM 80 L (3 CHARACTERS AT MOST) ASSIGNED BY USER Tn ThE T-TH 
NETWORK BRANCH IN THE DATA BRANCH list. THE ELEMENTS VALUE ♦ 

IF NOT SPrCIFIED IS REPRESENTED BY THIS SYMBOL 
TYPR (NtJMC) STYPX ( I > WMERE NUMC©NUMX(I) 

(USED ONLY FOR NETWORK TREE BRANCHES) 

TYPE (NUMC) “TYPX (LINK) WHERE NUMCsNUMX ( L I NK ) 

(USED ONLY FOR NETWORK BRANCHES NOT IN THE TREF) 

TYPX ( I ) aSPECIFl^S THE ELEMENT TYPE OF THE I-TH NETWORK BRANCH 

IN THE DATA BRANCH LISf« (MUST BE E ♦ I * 0 « R ♦ L » C » CC ♦ CV . VC» OR VV 
AND MUST BE COMPATIBLE W I TH THE UN I TS OF Th^ ELEMENTS VALUE) 
VAL (NUMC) sVALX ( I ) WHERE NUMCsNUMX(I) 

(USED ONLv FOR NETWORK TRfE BRANCHES) 

VALX(I) sELEMENT VALUE (E12.5) OF I-TH IslETWOHK BRANCH TN' THE DATA 
BRANCH LlsT 

wEIGT (I)®CONSTA nT TERM ASSOCIATED WITH THE VALUE OF THE SFG BRANCH T 


i 
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IV-2 . Flow Charts 

Program SNAP is divided into the following sections: 

Program MAIN (Subprograms 1 thru 12) 

Subroutine SFG (Subprograms A thru J) 

Subroutine FTREE 
Subroutine TREP 
Subroutine ARRAY 
Subroutine DECODE 

As indicated above, program MAIN is further broken down into 12 subprograms 
and subroutine SFG is divided into 10 subprograms. 
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Subprogram KAIH-1 

This program reads in some preliminary network data. 


Read in 

(a) problem name 

(b) HOD, NOB, KBASIS, LISTG, 
LISTC, LISTP , NIN, NOUT, 
NODA, NODB 

Set KBASIS to 8 if a zero valve 
has been read in. 

Write out the above information 
for reference purposes. 


* 

To Subprogram MAIN- 2 



Subprogram MAIN- 2 

This program generates the SFG routing matrix, creates a code for each 
symbol (excluding s), and sets up arrays for the constants and powers of 
$ associated with the branch weights. 
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SMBOL(IBO) - SYHBUL(LOB) 


307 

* 




Test 

NEST(LOB) 


KODE ( J , L) ■« 0 


SFG branch LOB 
contains a 
symbol 


SFG branch LOB 
does not contain 
a symbol 


AkL 


307 

* 


Test 
IG(L) 

where L » NLAST(LOB) 


^ K0DE(J,L) ° XG(L) 




A code has 
already been 
assigned to 
branches entering 
node L 




KP ■ KPH 
Initially KP=0 


&L 


Compare 
SMBOL(IBO) to 
iSMBOL(KP) 


not 

equal 


* 

952 


equal 




307 

A 


Compare 

/ KONSO(LOB) and KONSO(LOBX) 
\ where LOBX - MIX(KP) 


equal 




KODE (J ,L) « IG(LX) 
LX « NLAST(LOBX) 


not 

equal 


not equal 

the symbol associated 
with branch IBO is the 
inverse of the symbol 
associated with branch KP 


Compare 
KP to I BO-1 


752 


equal 
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2® 
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Subprogram MIN- 3 

This program codes the nodes of the SFG, and prepares the counters for finding 
all paths and/or circuits. 



To Subprogram MAIN-4 
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Subprograa MAIN-4 


This program finds all paths from node NIN to node NQUT and/or all circuits 
of the SFG. 


Preliminary 

(PF1-1) 

IFL0W(Il)-0 

Il-l,NNG 

K0NG(Il)-0 

Il»l,NNG 

NOP^KLAS 


KLAS^O 





I » 2 




jL 


Prelitninary (PF1-2) 

JX(1)«NIN 

J°H1H 

N?(1)°NIN 

IFLOH(NIN)°l 

IFL0W(K0UT)=>-1 




K “ KONC(J) ^ 
~ 


PF2(FIHD NEXT NODE) 
NP(I) n N(J,K) 

"ZZ 3ZZZ ZZ 

PF3 

(TEST ROUTING MATRIX) 
N ( J , K) 


< 0 


All loops thru 
a particular node 
found. Go to 
subprogram MAIN -5 to 
eliminate this 
node. 


100 

(Subprogram MAIN- 5) 


23 


25 


> 0 

continue to 
flower check 


* 

34 


All paths out of 
node J checked. 
Return to 
previous node 


25 

A- 


IFLCU(J)« b 0 
K0NC(J)-0 
J®NP(I-2) 
K0NC(J) s K0NC( J)+l 
I“l-1 
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another loop 
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Subprogram MAIN -5 


This program determines if circuits are to be found and if so modifies 


the SFG by eliminating node J. 
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Subprogram MAIN-6 
nd 

This program finds and stores all 2 order nontouching loops 
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increment 

LIR2 



End of 
Do loop 


NOONOOl 
NOTCH (HOC) “LI R2 

The following defines arrays NOTCH and NOCTOT: 
Consider the integer set 

Cl 3“ Cl ,2, . . . ,N2 } 

where N2 « number of nontouching pairs of loops 
Now consider the following subset of Cl 3 
Cl }“lN0CT0T(K>l)+l,N0CT0T(K-l)+2, . . . ,NOCTOT(K) } 


Then the set 


Cnotch(j) : js [i g 33 


is the set of loops that do not touch loop K 
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Subprogram MAIN- 7 

This program finds all nontouching loops of order greater than 2, and 
stores the associated code, power of 8, and constant term. 

The matrix ISET is given below In its general fora to aid in understanding 
the flow chart of subprogram MAIN-7. 

ISET(l,l)ISET(l,2)*«-ISET[l,NUP(l) ]***ISET[1,JAC(1) ] 

ISET(2, 1)ISET(2 , 2) • • *ISET[2 ,NUP(2) ]• • *ISET[2, JAC(2) ] 

• • 

• • 

• • 

ISET(KAP,1)ISET(KAP,2)“ •ISET[KAP,NUP(KAP) • *ISET [KAP, JAC(KAP) ] 
ISET(KAP#3,NUP(KAP+1) ]•» • 

where ISET(J,I), I-l f 2, . . . , JAC(J) 
is the subset of 

tlSET[j-l f mjP(J-l)+l],ISET[j-l,NUP(J-lH2] ISET [.J-l , JAC(J-l) ]} 

which does not touch the loop 
ISET[J-1,NUP(J-1) ] 
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All higher order 
loops not touching 
loop KIP have been 
found 




TEST 

KAP 


> 0 


490 


KAP-KAM ^ 
KUP(KAP)“0 p 


JAC(KAP+l)=-0 

HUP(KAP) a NUP<KAP)+l 


429 


* 
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CALL ARRAY (KAP+2 , TCOM SG , KXPOG , KSYMG , POLY , LI C , KIK) 
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increment 

MAPI 


*» Es\£] 


Update column counter of row KAP+l 
of ISET and insert the last loop 
found into ISET 

JAC (KAP+1 ) “J AC (KAP+1 )+l 
I SET [fCAP+1 , J AC ( KAP+ 1 ) ]=I SET (KAP , MAPI ) 


TEST 

JAC(KAP+l)-2 


^0 


425 

J 


< 0 


Row KAP+l of ISET 
has now been found 
and it has more than 
1 non-zero entries. 
Thus, procede to 
find the entires 
of row KAP+2. 


TEST \> 0 

J AC (KAP ) -NUP (KAP ) - 1 / 


429 

J 


440 


< 0 


Row KAP+l of ISET 
has now been found 
but it does not contain 
more than one non- zero 
element. Further, all 
loops of row KAP have 
been exhausted. Thus, 
it is necessary to 
back up one row and 
re-evaluated row KAP 


Row KAP+l of ISET 
has now been found 
but it does not 
contain more than 
one non-zero entry. 
Since not all loops 
of row KAP have been 
exhausted, increment 
KUP(KAP) by one and 
re-evaluate row KAP+l 
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490 

End of 
DO loop 
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Subprogram MAIN- 8 

This program decodes composite codes representing nontouching loops and 
sets up tags for use in printing out the symbolic transfer function. 


Set 


J3 ra l,NT0; J4-NSPTU 


P0LYU(J1 , J2)“0 Jl°l ,ITEXPS ; J2«NTO 

SEMPON ( J 3 , J4 ) a STA R( 1 ) ~\ 

SEMPOD ( J 3 , J4) ra STAR( 1 ) \ 

SIKB0N(J3, J4)“SB f 

SIMBOD(J3,J4)“SB 
NA(J5)°0 1 
NB(J5)°0 / J;> 

SB and STAR(l) are obtained from "data” statements 


L 

/ 

DO 646 JZ=1 ,LIL-1 
(LIL-1 different compo 
have been found in Sub 
— 

site symbol codes 
program MAIN-7,) 
— 


-k- 


KODY*»KSORT( JZ) 
ITOP(JZ) m O 


There is no syiabol 
associated with this 
value of JZ 



> 0 


y 

646 
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CALL DECODE (KOO,KBASIS,KODY,IZ,FB,JZ,SEMBOL,KQDF, 

KODI , ITOP ) 

This subroutine (a) sets 

A if terms having the code KSORT(JZ) 

ITOP(JZ)® I belong in numerator of output polynomial 
0 if terms belong in denominator 

(b) finds the set 
KODI(I), 1-1, 1Z 

where SEMBQL[K0DS(I) ], I®1,IZ are the corresponding 
set of symbols, and 

(c) finds the multiplicity of each individual symbol, 
SEMB0L[K0D1(I) ], and records these values in the 
array KODF(I), I°1,IZ 




DO 645 NZ=1, 1Z 


JL 


mQ 


The symbol corresponding 
! to KODI(NZ) is not to be 
j inverted in the output 


\ 

TEST \ «1 

kons[kodi(kz>] 


The symbol corresponding 
to KODI(NZ) is to appear 
inverted in the output, 
l.e. 


SEMBOlCkGM(NZ) ] 


NAK“NAK+1 

SIMBOH ( J Z , NAK ) ■ SEMBOL [KODI ( MZ ) ] 
SEMPON(JZ,NAK)°STAR[KODF(HZ) ] 


KAT“HAT+-1 

SIMBOD(JZ,KAT)°SEMBOL[KODI(NZ) ] 
SEHPOD(JZ,NAT)«STAR[KODF(HZ) ] 


NA(JZ)°NA(JZ)+1 


J NB(JZ)°HB(JZ)-H 
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Subprograa MAIN -9 


This program separates POLY Into the arrays POLYU and POLY for use 
in printing out the constant terras of the transfer function. 



Subprogram MAIN- 10 








-71- 


Subprogram MAIN-10 


This program normalizes the transfer function so as to have all 
positive powers of s 


522 


522 


522 


\ 

* 

DO 522 KAR°1, KIK-1 

i 

< 


> 0 


TEST 

MSORT(KAR) 


> 0 


< 0 


t- 


The power of s corresponding 
to row KAR of POLYU or FOLY 
is negative 


TEST 
MAXIMM-MSORT(KAR) 


< 0 




MAXIM is to be the 
absolute value of the 
most negative power of g 


MAXI Mf-MSORT (KAR) 


fEnd of 1 
LDO loopJ 


522 




■7 Continue 




^ increment 
KAR 


Using MAXIM make powers of 
s either positive or zero, 
i « e « 

MSORT(K) “MAXIM4-MS0RT (K ) 
K«1,2,...,KIK-1 


Subprogram MAIN- 11 
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Subprogram MAIN- 11 

Write out the matrix of constant coefficients for the numerator polynomial 
of the transfer function. Also write out the eyehole and e powers that 
correspond respectively to the columns and rows of the array of constants. 


Subprogram MAIN-12 

Write out the matrix of constant coefficients for the denominator 
polynomial of the transfer function. Also write out the symbols and s 
powers that correspond respectively to the columns and rows of the array of 


constants 



Subroutine SFG(NFI RST.NLAST.IXPON, WEIGHT, SYMBUL.EjQHSO, MIX, HEST, LIST, NIM, 
NOUT , KOD , KOB , LISTG, KODA , KODB) 

This subroutine generates a signal flow-graph (SFG) for the given network. 
The program io subdivided into subprograms A thru J. 
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Subprogram "A" 


This program uses DATA statements to define certain variables, nulls 
arrays, creates the SFG feedback branch, reads In network branch Information 
and call 8 FT REE to find a tree. 




Use DATA statement to define the following 
variables: Y,G,C,XQ,R,CL,Z,E,CI,CC,CV,W, 
VC , FB , ONE . 

\ 

/ 

Set NS(IC,IK)“0 \ 

NF(IC, IK) 53 0 J ^ A * 
NEST(IG)°0 \ 
KONSO(IG)“OJ ’ 

INTREE (11 )*»0\ , 

JROW(Il)oO J 11-1 » 
MCM) , L0=0 , LI ST«*1 , LINK^Q 

1 

NNG; XK°1,NNG 

NBG 

NNG 

jy 

U. 


Create SFG feedback branch used to make 
the SFG "closed". 

NLAST(1) e NIN 

IXPON(1)«*0 

WEIGT(l) a -l 

SYMBUL(1)«FB 

KONSO(l) a O 

NEST<1)=»1 

Note: NFIRST(l) Is determined In Subprogram "I" 
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V 


Subprogram "B" 
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Subprogram "B" 


This program sets up tree branch information and creates a routing 


matrix and sign matrix for the tree. 

i 

DO 21 NU°1, NOD-1 


& 

Put network tree branch information in terms 
of the branch labels specified by user. 

Let IO^INTRE(HU) »NUMC°HUMX(IO) 

Then TYPB(HUMC)“TYPX(XO) 

JB(NUMC)°JBX(XO) 

LB(NUMC)°LBX(XO) 

S YI1 ( HUKC ) D S YI-1X ( 10 ) 

XQUAL (NUMC) **1QUALX ( 10 ) 
VAL(NUMC) a >VALX(I0) 

HUML(HUMC) “HimLX(XO) 

INTREE (NUMC)*»1 



NE. 
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Increment 


IB(JF,LF) E HUMC 
IB(LF, JF) a NUMC 


The following generates a routing 
matrix, HP, and a sign matrix, 

HS, from the chosen network tree 
for use in finding the SFG 


±. 


JF=JB(NUKC) 
LF ra LB(HUMC) 
JROW(JF) ra JKOW(JF)+l 
HF[JF,JEDVJ(JF) ]= LF 
NS[jF,LF>l 
JR0 W(LF) ei JR0W(LF)-}-1 
NF[LF,JR0!4(LF) >jf 
ns[lf,jf>-i 


'End of " 
21 LdO loop. 



. — 1 

Continue 


\ 

f 


Subprogram "C" 
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Subprogram "C" 


This program generates SPG information from branch node to link node. 


15 1 
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Subprogram "D" 


This program generates SFG Information from link nod© to branch node. 










149 

Subprogram "C" 
* 


KFIRST(LIST)= s XFXM 
NLAST (LIST) “INIT 
SYMBUL(LIST) 0 SYM(INIT) 


NPLA=NPL“ 1 -LON 

N 

' 


TEST 


0 


151 

Subprogram "C" 
* 


NPLA 
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Subprograa "E" 

This program sets up SFG information for VC type control sources. 


165 NUMB was found from 







- 85 - 




branch NUMB 
is impedance 
type 


COMPARE 

IQUAL(NUMB) to 

IQ 


TEST 

KUNO 


NEST(LIST)-! 
WEIGHT(LIST)°1 . 


Network branch 
NUMB contains 
a symbol 


WEIGHT(LIST) 
■ CVAL(NUMB) 


HEIGHT(LIST)“ /CVAL(NUMB) 


Using the newly created 
SFG node representing 
the controlling voltage, 
a SFG node can now be 
generated for the dependent 
current source 
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123 

Subprogram "C" 
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Subprogram "F" 

This program sets up SF6 information for CC type control sources. 


265 

4 

The structure of subprogram "F" is 
completely analogous to subprogram "E" . 


V 

123 Subprogram "C" 


Subprogram "G" 


This program sets up SFG information for VV type control sources. 


360 

i 

This program is cycled over all voltage 
controlled voltage sources in the network. 
A list of these sources is maintained in 
the array 

IW(ffl) , MI-1, MO 

Other than the above, the structure of 
subprogram "G" is completely analogous to 
subprogram "E" . 




460 Subprogram "H" 
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Subprogram "H" 


This program sets up SFG information for CV type control sources. 


This program is cycled over all 
current controlled voltage sources 
in the network. A list of these 
sources is maintained in the array 
ICV(MI),M3>1,L0 

Other than the above, the structure 
of subprogram "H" is completely 
analogous to subprogram "E" 


V 

515 Subprogram "I" 

The SFG is now 
complete except 
for the output 
node. 
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Subprogram "I" 


This program generates the output node of the SFG. 
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I 


514 



Subprograa "J" 
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Subprogram "J" 

This program orders SF6 Information for input to main program. That 
is, a mapping function MIX is found which reorders the SFG information so 
that the routing matrix N(J,K) as calculated in MAIN will automatically 
have its entries decrease as K increases. 










KIX(IU) and 
j MIX(IL) Exist be 

I interchanged 


MXL ra MIL(IL) 

MX(IL)°MIX(IU) 

MIX(IU) n MXL 
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Subroutine FTSEE(TYPX,JBX, 
LBX , INTRE , NOTREE , NOD , HOB) 


This prograa finds a tree of the network to use in generating a SFG. 









K“Mf 1 

initially M“0 


branch 
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\ 


Subroutine TREP(NIN,NOUT,NF,NP,NPL) 


Given the routing matrix for the network tree with input and output nodes 
specified, this subroutine finds a node list representing the path between 
these two nodes. 


25 

t 


J, 


Preliminary 


Set JX(I5)~0 
NP(I5)>=0 
JMEM(I5)=-0 
KMEM(I5)=0 


20 

25 


15-1, HNG 


/ NPL-0, JX(l) ra JX(2) ra NIN, 
I"l.tT B, NIN,NP(l)-NIN 


1 

K* 

*0 


t : 

) K^K+1 



TEST 

NF(J,K)-N0UT 


["Path ”| 
Lfinished J 


“0 





4 

50 


/ TEST \ “0 

\NF(J,K) / 


TEST 

J-NIM 


= 0 


> 0 


[Flower check] 


js k 


[Flower formed] 


TEST 

NF(J,K)-JX(I) 


* 0 


V 

60 

[Backs tep ] 


100 

Stop 
no path 
found 


*0 




45 
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20 

11 


45 

/ 


Store and resaember vertex 

NP<1)«NF(J,K) 

JMEH(I) a J 
JX(I+1)“HF(J ,K) 

J=NF(J,K) 

KMEM(I) e K 


25 

Backs tep 



J“JMEM(I) 



K°KMEM(I) 

\ 

■ | 

I°I-1 



Final path vertex and 
path length 
NP(I+1)=N0UT 
NPL«I+1 


k- 


50 


Return 




100 



- 97 - 


Subroutine ARRAY ( J SI G , XCON , J XPO , JKOD , POLY , LIL , KIK) 

This subroutine takes the constant associated with each loop or 
nontouching combination of loops in the SFG and stores it in the matrix 
POLY. It does this by comparing the code and exponent of the given loop 
combination with the codes and exponents assigned to the columns and rows 



10 j* 0 
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V 

20 


JKOD-KSORT(NH) 7 



! KSORT(LIL)“JKOD 

24 y NNX°LIL 

! LIL^LILrfl 





Return 
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Subroutine DECODE (KOO , KODY , IZ , PB , JZ , SEMBOL , KODF , KODI , ITOP , KBASI S ) 


This progress decodes the composite symbol codes. 

k 


IZ*»0 

1MCBASIS-1 


ii— 


DO 3 J«l,KOO 





I POWE R® M . AND . KO DY 


1 


TEST 
I POWER 


£0 


M. 


KODY contains 
code for 
SEmOL(J) 


COMPARE 
SEMBOL(J) to FB 


1 


.NE. 


.EQ. 

The tern having 
code KODY belongs 
in numerator of 
transfer function 



IZ»IZ+1 

j K0DF(IZ)«=IP0WER 
K0DI(IZ)®J 

t 

1 - A- 


3 pnd of " 
LDO loop. 

/ 

increment 

J \ 

k 

t 

KODY®KODY/KBASIS 






IT0P(JZ) a l 





n o o o o o o 
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IV-3. Program Listing 


SNAP 

This PROGRAM FImUS THF SYMbOLlC TRANSFER 
FUNCTION Ok IMMITANCE FUNC1JOM OF A 
LUMPED LlNfcAR TIME. TNVARIAinT NtTwOPK 




C The. FOLLOWING arrays ARE ASSOCIATED WITH the nftwork 
C CHARACTERISTIC NhN U'EFINED in program MAlN-i) 

DIMENSION Lt ( 35) ♦ IG < 35) ♦SMROL (35) 

DIMENSION Ir Low ( 35) ,NP (3s) ♦KOnES (35) .KONC (35) 

DIMENSION N*( 35 ♦ 3b ) * CONS ( $5,35) ♦ RODE ( 3.5 ♦ 3b ) »IaP0(35»35) 
G*****#*»**tt**«****« **»*#*•»*#*******•»******«****»#**#*******##*****#*** 


C The FOLLOWING aHrAYS ARF ASSOCIATtU WITH THE NF T WORK 
C CHARACTER I ST IC HbG 

DIMENSION NfIRSI ( 100) ♦NLaST ( 1 0 o ) t IXPON ( 100) ♦wElGT ( 100) 
DIMENSION SYMBOL (lOQ)'Mlx(100)t OVAL (100) 

DIMENSION RqNSOUOO) ♦ NEST (100) ♦TYPE ( 1 0 0 ) 


c thf following arrays are associated with the mfiwork 

C CHARACTF RIST IC NPAC 

UlMENSlON CoNST ( 300 ) ,K0 DeT ( 30 0 > ♦TXPOT ( 300 > ,MaPO (300 ) 
DIMENSION inqCTOT (3C0) ♦NUk ( 300 ) ♦ JaC ( 300 ) 

DIMENSION NdCOOE ( 30 0 ) 


C Thf FOLLOWING ARRAYS ARE ASSOCIATED WITH THE NEIwoRK 
C CHARACTERlST IC s nTo* NSP T ♦ AND NpXPS 
DIMENSION N a ( 1 b 0 ) ♦ N B ( 1 5 0 ) 

DIMENSION t\ 0 NS (FU) i K O D I (<?C) ♦SfMbOL (?0) ♦ KOUF ( £ 0 ) 
DIMENSION MSORT (15) ,RSORl (150) ♦POL YU ( 15, IbO) 
DIMENSION PqLY ( lb, 150) ♦ I TOP , 150) 

DIMENSION SIMMON ( 1 5» 0 ♦ 1 0 ) , S I MBnD ( 1 50 ♦ 1 O ) 

DIMENSION SEMROn ( 1 50 ♦ 1 0 ) ♦ SEMPOU (1 50 , 1 0 > 



?( h 

it I 

50 

7 
HQ 


0 


2 3 0 


Q 


3 c f 
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C***«HMK»*»*«*****«*#*«****«********************<HHHHMHHHHKHHHHHMMHHHHM»« 390 

<>0 0 

C The FOLLOWING a&RAYS ARE ASSOCIATED WITH THE NF I WORK 410 

C CHARACTERISTIC*; nRI» NCI» NEON ♦ AND NKS 420 

DIMENSION UET (1S»100) *NoTChU200) »STAR(9) 430 

COMMON SEMPnN*SEMPOD 9 PULY 450 

COMMON/CI/MsORT»KSORT 460 

COMmON/C2/NNG*NHG 470 

COMMON/C3/NtXPS»NTO 480 

C0MM0N/C4/NSPT 490 

EQUIVALENCE (IXHOdfl) »NOTCH(l> *SIMBOKI(l«l)) 500 

EQUIVALENCE (CON$d«l) tlsET (1*1) «SIMB0D(U1) ) 510 

EQUIVALENCE (KOUEdtl) •PoLYU(ltl) ) 52o 

DATA 0 A$h/2m // 530 

DATA Fti « SH/3H FH*3H \ / 540 

DATA STAR(l) »STAH(2 ) *STAh( 3)/3H * 3H*«2 t 3 h**3/ , 550 

DATA STAR(4) »STAR(5) ♦ S T Aw ( 6) /3H<**4 » 3H**5 » 3H**6/ ' 560 

DATA STaR(7> *STAK( 8 ) ♦STAn(9)/3H*#7»3H#*8»3h*«9/ - 570 

DATA ONE/3H 1/ 580 

590 

600 

C PROGRAM MAIN— I 610 

C PRELIMINARY INPUT INFORMATION 620 

630 

640 

C NRNrNUMBER OF BRANCHES IN NETWORK 650 

C nBGbNUmBER OF BRANCHES OF SFG 660 

c NTOsNUMBER of Terms in output 670 

C NSPTsNUMHER Op SYMBOLS PER TERM 680 

C NExPSsNUMBEH nF DIFFERENT POWERS OF S 690 

C NPACSNUMHER or PATHS PLUS CIRCUITS 700 

C nrt=maximum number OF nunTouching LOOPS 710 

C NCI=MAXIMUM NUMBER OF LOOPS not touching any given loop 7?0 

C NF:0NSNUMB£R Op nOnTOUCHING PAIRS OF LOOPS 730 

c NRS=NUMAER OF REPEATED SYMhOLS (NUMBER OF NETWORK , 740 

C ELEMENTS ASSIGNED SAME SYMBOL) > 750 

nBN*35 760 

NRGslOO 770 

NTO= 1 50 780 

NSPT=20 790 

NF XPS= 1 5 800 

NPAC-300 - 81f> 

NR 1*15 ' , 820 

NC 1 = 10 0 '* '•.< 830 

NEON* 1200 840 

NPS*9 850 

c NSPTU=NUMBER nF SYMBOLS IN NUMERATOR OF EACH TE«M 860 

C NRTGaNUMBER UF BRANCHES IN TREE UF SFG 870 

C NNG=NUMHEW OF NODES IN SFG 88 q 

NNGsNBN 89o 

NSPTU=NSPT/2 900 

N8T G=NHN 910 

CALL StCONU (TCOMP) 920 

WRITE (6 » 160A) TCOmP 930 

1600 FORMAT (1X»48h CUMPILATIqN T I Mt IN SECONDS , F 1 5 . 8/ ) 940 

1111 CONTINUE " 95o 

WRITE ( 6 * S 1 9 ) 960 

519 FORMATllHl) 970 

CALL SECOND (TSTART) 98o 

C THE NEXT 6 CARDS ARE FDR PROBLEM TDEnT IF IC«T ION ON THE 1ST DATA CARD 990 
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RE AD (5* 1150) (WEIGT(J) *J= 1*72) 

1150 FORMAT ( 72A i ) 
lF(E0Ff5)lllll»11112 

11111 STOP 
nil? CONTINUE 

WRITE (6tll60) (WEIGT(J) * J=l*7l> 

1 ISO FORMAT ( 1X,7 1AI//) 

00 1151 Jsl ’ 7;f 

1151 WE I6T ( J) =0 • 

RE AD ( 5 9 1240) nOD*NOB*KBASIS*LISTg*LISTC*LISTP 
1240 FORMAT (315*5X931 1) 

IF(KHASIS) 1357* 1357*1358 

1357 KRASIS=8 

1358 CONTINUE 

RE AD ( 5 9 1 ) NiN*N0UT*n0DA»n0D8 
1 FORMAT (415) 

WRITE (6, 720) MOO 

720 FORMAT (2X* loNNUMHER OF N0DE8=,13l 
WRITE (6*72l)NQ8 

721 FORMAT (2X» lqHNUMbER OF Bh ANCHES= *13) 

IF (LISI 0) 72^» 723*722 

722 CONTINUE 
C LIST SFG 

723 IE(LISTC) 725*725*724 

724 CONTINUE 

C LIST ALL CIRCUITS 

725 IF (LIST?) 72if * 726*727 

727 WRITE(b*72H) 

728 FORMAT (2X*24NLIST ALL PATHS FROM NOOK *I3*2X* 7HTO NODE *13) 

726 WRITE (6*729) NlN 

729 FORMAT (2x,25HELEMENT NUMBER OF SOURCf=*I3) 


IF (NOUT ) 1B0?». I 802* 1804 
1804 CONTINUE 

wRITF. (6*730) NOUl 

730 FORMAT (?X*3«hELc.mEnT NUMBER 
GO TO 1 80b 

1 802 CONT INUE 

WRITE (6*731) NOUA 

731 FORMAT (2X*33HR0S1TIVE OUIPUT 
WRITE (6*732) NOUB 

732 FORMAT (2X.33HITEGATIVE OUTPUT 
1806 CONTINUE 

WRITE (6,850) KHAS1S 
H50 FORMAT (2X*22HriASE FOR SYMBOL 


ASSOCIATEO wITH 0UTpUT=*I3) 

VOLTAGE TERMINAL'* 13) 
VOLTAGE TERMINAL'* 13) 


CODES'* 14) 


program main — 2 

T A*E SFG BRANCH INFORMATION AS FOUND 
BY SUBROUTINE SFG ANU GENERATE 

(1) routing matrix information 

N( J*l\> * AND LI ( J) 

(2) SFG BRANCH VALUES IXP0(J,L)* CONS(J*L), 

KOUE(J*L) WHERE JsNFIRST(T)* L=NLAST(I), AND 
IsBRANCH NUMBER 

TOGtlHtR WITH THE SYMfciUL StMBO(.(K)* K=l,?* e ..*Ml 


call subroutinf to formulate the signal flow gpmph* sfg 

Call sfg (NF1RSI*NUASI * I XPON* wE 1 G T ♦ SYMBUL * KOnSO * M 1 X , NEST , | 1ST* 
ININ* NOUT *NOU»NOB»LlSTG*N(jOA*NnOB) 


i 0 4 rl ] 

10 50^ 
1060 
!07rH 
lOK.^J 
1 09n 

1 1 0 
11 K 


1130 


1370 


15? l 

1 5 Sri 

lbyn 

16 on 
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IF (NOB) 1111*1111*1920 1610 

1920 CONTINUE 1620 

CALL SECOND (T2) 1630 

WRITE (6*93;) 1640 

TSFG=T2-TSTflRT 1&50 

WRITE (6,16o2) TSFG 1660 

1602 FORMAT (1X,2(HTIM£ FOR FORMULATING SIGNAL/ 1670 

122H FLOW GRaPH IN SECONDS * F 15 . 3/ > 1680 

IROsO 1690 

K0=0 1700 

MICH=1 1710 

K=0 1720 

MG=1 1730 

JL AS* l 1740 

NCIRsl 1750 

ININ=NIN 1760 

INOUT*NOUl 1770 

KOOsO 1780 

DO 301 iNKsV *NSPT 1790 

301 KONS(INK)=0 1800 

DO 300 INK=l»NNG 1810 

300 IG ( INK ) =0 1820 

1830 

C FIND IXP0(J,L) »CONS(J«L> 1840 

GO TO 307 1 850 

305 MG=K8AS1S*MG i860 

MICH=MICH*1 1870 

307 I BO- 1 80* 1 1880 

IF (LIST-IHO) l9»4»4 1890 

4 CONTINUE 1900 

L08 = M1X< [HO) ' ~ ' 19 VO 

JsNFIRST (LOB) 1920 

L=NLAST (LOG) / 1930 

IXPO( J»L) =XxPON(L08) 1940 

CONS ( J*L> =WfIGT (LOB) 1950 

I960 

c. FIND ROUTING MATRIX - .... 197n 

8 IF (J.EQoJLAc;) GU 10 10 i9Rn 

LT(JLAS)=K 1990 

K1=K*1 , 2000 

IF (JLAS-NIN) 2R»27»28 2010 

27 N ( JL AS » K 1 ) =- 1 20?o 

GO TO 29 2030 

28 N(JLAS»K1)=0 2040 

29 JLAS= JL AS* 1 2050 

K=0 2060 

GO TO b 2070 

10 K=K*1 2080 

N(J»K)=L 2090 

2100 

C FIND KOUE ( J » L ) AND SEMBOL(KUO) 2 1 1 0 

SMH0L( lBO)=c:YMriUL(LOB) 2120 

mOOE = NEST (Lob) ' 2130 

IF (MODE) 335»3lfa»335 2140 

335 IF (I^(L) )5»960»b 2150 

5 KODE ( J ♦ L) = lb ( L ) 2160 

GO TO 307 2170 

960 CONTINUE 2180 

KPU= I BO- 1 2190 

IE (KPU) 953*gb3» J15 2200 

315 00 952 KP=1,KPU 2210 
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IF(SMBOL(IBO> .NE,SMROL (KP) ) GO TO 952 
LORX = MU (KP) 

IF (KONSO(LOr) -KUNSO(LO bX) ) 952,956,952 
956 LXsMI.AST (LOfl*) 

NODE ( J»L) *1r (LX) 

GO TO 307 

952 CONTINUE 

IF(SMBOLUbo) .EO.ONE)GO TO 3U 

953 IG(L)=MG 
KOO=KOO+l 

SFMtiOL (KOO) 2 SMdUL ( I BO) 

KOOE<J*t)=Ic(L) 

IF (KONSO(LOr) ) 3,3,2 

2 KONS(KOO)=l 

3 CONTINUE 
GO TO 305 

316 KODE ( J * L ) =0 
GO TO 307 
19 LT(JLAS)=K 
MUK+l 
N ( JLAS * K1 1 ) = U 


PROGRAM MALN--3 

NULL CtHlAt.N ARRAYS, SET COUNTERS, AND DEFTNE 
A CODE FOR EACH NODE OF THE SFU 


MPL = 0 
KTKsl 
L T L s 1 

DO 601 KAM = 'i,NEXPS 
DO 601 KIM=1 »NT0 

601 POLY (KAM,KlM) =0 

DO bO? K P 1 = 1 » N E X P S , ^ 

602 MSORT(KP1)=u 

00 950 K02- 1 » NSP T 

950 KOD I (K02) =0 

DO 603 KP2 = 1 * NTO 

603 KSOR T ( KP2 ) s n 

IR = 1 

NF 1 R= 1 
KN0 = 0 

KOOES( 1) =1 
00 2000 JS=2*NNG 
2000 K0(>ES<JS)=2*NQUtS< JS-1) 

IF (LISTP) 175*175,1116 
1116 WRITE (6, 170) NINtNOUT 

170 FORMAT ( 1 7H 6ATmS FROM NOuE 12, 9H TO NODE 12//) 
WRITE (6,190s) 

1905 FORMAT (5X, l7 M NO. NODE LIST) 

175 CONTINUE ... 

IF(LISTP)11 1 3,U13.2J 
1 1 13 K3=LT (WIN) 

N (NIN»K3) =0 
K2=LT ( 1 ) * 1 

N (1 ,K2) s-l 
N I Ns 1 
NOUT - 1 

klas=o 

24 NFIR=0 


222(0 

223(U 

2240 

225rpj 

226 cfj 

227 0 
2280^ 

230'U 

2310 

232fR 

233 ( 

234 0 U 
235n 
236([1 
237<U 

235 0 
23P0-| 
240 f I 
2410* 

2420 
243ffl 
244fLj 
2450 
246 (V-j 

24 7 0 
245(P 
2490 
25*0fpj 

25 hu 

2520 

253 (jR 

254 0. | 

255 fH 
2560 
257 (p 

25Bf[J 

2590 
26()0_ 
26i 
2b2(H 
2630 
264 (H 
265rM 
2660 
2670-, 

268 f 
269(M 

2700 
2 7 1 rp 

2730 

275(0 
276M 
2770 
2 7 « .r? 

279([j 

2800 

2810 
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IF (LISTC) 12o9»l209*l219 2830 

1219 CONTINUE 2840 

*HITE(6,177) ' 2850 

177 FOPMAT (1X 9 8MCIHCU1TS//) 28(50 

WRITE ( 6 * 1905) 2870 

KNO=0 2880 

1209 CONTINUE 2690 

2900 
2910 

PROGRAM MA1N--4 - - 2920 

PaTH-rINDInG ALGORITHM 2910 

IN aODITIO'M, STEP PF7 CALCULATES THE COMPOSITE 2940 

COOE ♦ CONS f ANT « A NO EXPONENT OF THE PATH 29S0 

2960 
2970 

PFl (PRELIMINARY) 2980 

OO 1112 I ZOa 1 * NNG 2990 

1112 I FLOW ( I ZO ) =0 3000 

UO 31 1 1*1 »nng - - 3010 

31 KONC ( I l ) = 1 ' 3020 

nOP=KLAS 3030 

KLAS=0 ‘ / 3040 

23 1=2 y 3050 

J=NIN 3060 

NP(1)=NIN 3070 

IFLOW (NIN ) =1 3080 

1FL0W (NOUT ) =-l 3090 

3100 

25 K=K0NC(J) 3110 

3120 

PF? (FIND NEXT NODE) 3lJ0 

NP(I)=N(J,K) 31 VI 

3150 

PF 3 (TEST ROUTING MATRIX) 3i60 

I F ( N ( J » K ) ) 100*60*34 3170 

3 1 80 

PF 4 (TEST FOR FLOWER) 3190 

34 njksm(j,k) ; 3200 

IF ( IFLOw (NJ K ) ) 50*38,26 3210 

26 KONC ( J) =KONv- ( J) *\ < 3220 

GO TO 25 3230 

3240 

PF5 (PREPARE FoP NEXT NODE) 3250 

38 J=Np(I) 3260 

IFLOw (J)=l 3270 

1=1+1 3280 

GO TO 25 3290 

3300 

PF 6 (BACKSTEP) 33lo 

60 IFLOW (J) =0 33?0 

KONC ( J ) = 1 3330 

J = NP( 1-2) ; 3340 

KONC (J) =KON(~ (J) +1 3350 

1=1-1 3360 

GO TO 25 3370 

3380 

PF7 (FINISH PA T H) 3390 

50 KONC ( J) =KUNc ( J) +1 3400 

KLAS=KLAS+1 3410 

X 3420 

C Finn CODE FOR NODE PATH 3430 

/ 
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NPCODE (1R) 

1SU=I-1 

DO 2002 IS=1»ISU 
h(0DS = NP(IS) 

?0 02 NPCOOF. ( IR) =nHC0DE ( IR) ♦KOUES (NOUS) 

CALL ARRAY aNu WHl IF 

I F ( Nc I R* EQ» l ) WO TO 179 
IF (LISTC) 12 0 H»120B,1206 
1206 CONTINUE 
KRU=I 

179 KNO=KNU*l 

WRITE (6*110) KNO » (Mp (KK) ,KR = 1 ,t\RU) 

110 FORMAT (4X»i3*6X05I3) 

120ft CONTINUE 

IF ( NF I R . EU» 1 ) GO TO 320 

KODET(IK)=0 

CONST(lR)=l. 

IXPnT(lR)=0 

IENOsI / 

00 319 KEW=Z*IENO 

JNODE=NP(KEw-l) 

LNOOEsNp (KE w> 

KOOET ( I K ) sKOUET ( IR) ♦ KODE ( JNODF ♦ LNOOF ) 

CONST ( IR) =CoNST ( IR) *CONS ( JnODF ♦ LNODF ) 

IXPOT { IR) sIXPOT ( IR ) ♦ I XPO < JNODF. * LNOOE) 

319 CONTINUE 
CONEW=CONST ( IR) 

IXNEW=IXPOT ( IR) 

KONEWsKODET \ IR) 

CALL ARRAY ( 1 * CONEW ♦ I XNEd , KONEW % POLY ♦ L I L ♦ KlK ) 

320 CONTINUE 


IR=IR+1 

IF ( IR~NPAC) 1361 ♦ 136l»136t) 

1360 wRITE(6 ? 1362) 

1 362 F0RMAT(1X»39HN0. OF CIRCUITS EXCEEDS L 1 M I T- 1 (mCRE ASE » 
126HDINEnSIOnS CONTAINING NPAC) 

1361 CONTINUE 
GO TO 25 


PKOGRAM MA1N--5 

MODIFY THE SF^ BY REMOVING EVF.ftY BRANCH CONNECTED 
TO THE NOUt THROUGH WHICH ALL CIRCUITS HaVF JUST 
BtEN FOUND 

I 

100 T 3 = 0 » 

IF(NCIR-I)2ol0»l02»20l0 

102 CONTINUE 

IF (NFIH-1) 10^*2010»104 

103 K4 = LT ( N 1 N ) + 1 
N (NIN»K4> =0 
K 5 =LT ( 1 ) +1 

N ( 1 » K5 ) =•* 1 
N I Ns 1 

NOUT=l 
60 TO 24 


344 d. 

3460 
34 7 t 
34ft. f j 
349 (j 


JHho 

38<) A 

39 ( \ 'I 
39lli 

39?n 

39 

39STT 

396n 

Xfl 

3990 
4 0 M 
4 0 f I 
4 0 H 

40 10 
4 0 a 
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104 IF (NIN-JLAS) 1 05*200,200 - - 4050 

105 NT N = J* 1 4060 

NOUT=J*l 4070 

KONC ( J) =1 40H0 

NY = LT ( J ) ♦ 1 4090 

N(J*MY)=0 4100 

00 109 JCsNtN,JLAS 4ll0 

LCOL=LT(JC> 4120 

IF (LCOL.E.U. r\) GO TO 109 4130 

IF (N(JC«LCOL)-J) 109(107*109 , - 4l40 

107 N(JC*LCOU=0 . 4150 

LT ( JC) =LT ( JC> -1 4160 

109 CONTINUE 4170 

NZ»LT(NIN)+) 4160 

N(NIN*NZ)=-1 4190 

N()UT = NI N 4200 

GO TO 23 ’ 4210 

2010 CALL SECOND ( T 3 ) v 4220 

TPATHsT3-T2 - 4230 

WRITE (6*201-?) NOPtTPATH 4240 

2012 FORMAT { 1X» l^HTlME FOP FINDING, Il0»l7H PATHS»IN SFCONDS , F 1 5 . 3/ ) 4250 

IF {NCIR-I ) 250, 103,250 4260 

200 CONTINUE 4270 

NOL = KL AS 42«0 

CALL SECONU'(T4) 4290 

IF (T3) 2014,2020,2014 4300 

2014 TCIRsT4-T3 4310 

GO TO 201b 4320 

2020 TCI«=T4-T2 4330 

2016 WRITE(6,160^)NOL»TCIR ‘ 4340 

1 603 FORMAT ( IX, IoHTImE FOR F 1 ND I NG , 1 1 0 , 1 8H F IrSt Oki)EP LOOPS , / 4 350 

HIM IN SECOnDS«F15.3/) 4360 

4370 
4 3 GO 

PROGRAM MAIN— 6 4390 

FIND StCONu ORDER LOOPS 4400 

44lp 

4420 

NOL=KLAS 4430 

KHOL=0 4440 

DO 257 KOW=i>NPAC 4450 

257 NOCTOT (KOw) 70 446o 

LOW l=NOP * l ' - 4470 

NOC = 0 4450 

NOLl=NOL-l 4490 

DO 203 L I R 1 sLOW 1 » NOL 1 4500 

L0W2=L1R1*1 4510 

DO 202 LIR2=L0W2,N0L 4b?o 

NANsNPCDDEtLiRl > *AND.NPCUDE (LIR2) 4530 

IF(NAN) 202»20l*202 4540 

201 CONTINUE 4550 

TCONS2=CONST ( L I R 1 > *CONST (LIR21 4660 

KXP02=IXPoT (LIR1) ♦IXPOT (L1R2) 4570 

kSYM 2=NUOET (LIR1) ♦XOOET (LlR2> 45 «o 

CALL ARRAY (2» TCUNS2,KXP02*KSYM2,P0LY»LlL»KIK) 459 o 

KHOL = KHOL-fl 4600 

NOC = N(jC*] 4610 

IF (NOC-NEUN) 139b, 1396, 1395 46?o 

1 395 WRITE ( b , 1 39 n ‘ 4b30 

1397 FORMAT ( 1X,46HINCKEASE NEoN-THE DIMENSION OF THE ARRAY NOTCH) 4640 

1396 CONTINUE 4bSc 
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NOTCH(NOC)=Ul«H 

202 CONTINUE 

203 NOCTOT (LIH1 ) = NOL 
NOCTOT (NOL) =noc 

PROGRAM HA I N — 7 

FIND ALL LUOPS Of ORDER 

6HEATER THAN 2 

GENERATE ThE FIRST ROW OF ISF.T 
NIPL=N0F*1 
K, APMO X = 1 
INKOal 

DO 1170 ISCrNIPLtNOL 
INKl=NOCTOT ( ISC) 

IF ( ISC-1) Il7l»im»ll72 

1172 I NK 2 s NOCTOT ( ISC-1) ♦ 1 - ' >■ • 

GO TO 1173 

1171 1NK2=1 

1173 IF ( INKl-INKp-INKO) 1 170 » 1 i 70 » 1 1 7S 
1175 INK0=1NK1-INK2 

1170 CONTINUE 

IF(INKO-NCI) 1391 ♦ 1391,1 3R0 

1390 WRITE (6, 139?) 1NKO 

1392 FORMAT ( IX, 52H1NCREASE NCl-THE NO. OF COLOMNS IN DIMENSION OF TSFT) 

1391 CONTINUE 

DO 490 NIP=n1PL,nOL 

INKUsNOCTOT (NIP) 

1FINIP" 1)210, 210,211 

211 INKLsNUCTOT (NIP-1) +1 
GO TO 212 

210 INKL«1 

212 CONTINUE 

IF(INKU-INK| )490, 490,410 
4l0 JIP=0 

DO AH 0 INK=1NKL, 1NKU - . 

JIP=JIP*1 

460 ISET(1»JIP)=N0TCH(1NK) . “ 

MAPO (NIP) al N KU-lNKL*l 

INI ATE PROCESS FUR FINDING ' 

HIGHER ORDER LOOPS 
DO 430 K A I = i , NP AC 
JAC ( K A T ) sO 
430 NUP(KAT)=Q 

JAC (1 ) = MAPG (NIP) 

KAP = 2 

440 KAPsKAP-1 

IF (KAP) A90,490, A29 
4?5 KAP=KAP«1 

IF (KAP-NRI) 1350, 13S0, 1351 

1351 WRITE (6, 135?) 

1352 FORMAT ( IX, A9HINCHEASE NR 1 -THE NO. OF ROWS IN DIMFNSJON OF TSFT) 
1350 CONTINUE 

NUP (KAP) =0 
429 K AP 1 =K Ap ♦ 1 
JAC ( KAP 1 ) =0 
NUP (KAP) BNUp (KAP) +1 


512 


52 ' 
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C LABEL LOOP OF FIRST CM 5270 

NAP=NUP<KaP) 5280 

IF (KAPMAX-KaP) 1347 * 1348, J 348 5200 

] 347 KAPM&X=KAP 5300 

1348 CONTINUE 5310 

1SAT=ISET (KAPtNAP) 5320 

C ' ' 5330 

C TEST LOOP OF REMAINING CKTS b340 

MAPU = JAC (KAp) 5350 

MAP|. = NUP(KAp) *1 '• 5360 

00 435 MAPI=MAPL»MAPU b370 

ISOTsISET (KaP*MAPI ) 5380 

KAN=NPCOOEUSAT) .ANO.NPCOOE ( ISOT) b39 0 

IF (KAN) 435*455,435 5400 

455 CONTINUE 5410 

C 5420 

C WRITE 5430 

TCONSG-CONST (NIP) 5440 

KXPOG=IXPOT (NIP) 5450 

KSYMGsKOOET (NIP) ' 5460 

DO 477 LPO=l»KAP 5470 

ITIC=NUP (LPO) 5480 

ITUCHsISEr(LPO»ITIC) 5490 

TC0NS6sTCUN(;G*CUNST ( ITUCm) 5500 

KXPOG=KXPOG*1 XPOT ( ITUCH) ' 5SlO 

477 KSYMG=KSYMtt*K00t 1 ( ITUCH) 5520 

TCONSG=TCONsG*COnST ( ISOT) 5530 

KXPOGsKXPOG*I XPOT ( ISOT) 5540 

KSYMG=NSYMG*KOOET ( ISM ) 5550 

KAPPrKAp*2 5560 

CALL ARRAY (kAPP»TCONSG»KxPOG»kSYMG»POLY,LIL»KIK) 5570 

KHOLsKHOL* 1 5580 

C SET COUNTERS 5590 

423 KAPUKAP+1 5600 

JAC(KAP1)=JAC(KAP1) +1 5610 

JACK=JAC (KApl) 5620 

ISET (KAP1 , JmCK) =1SET (KAP,MAPI ) ' 5630 

435 CONT INUE 5640 

JACK=JAC(kAp1) 5650 

IF(JACK-2)4il»42b*425 5660 

431 IF (JAC (KAP) -NUP (KAP)-l) 440*440*429 5670 

490 CONTINUE 5680 

CALL ARRAY (2* 1. *0*0* POLY* LIL*KlK) 5690 

CALL SECOND (I 5) 5700 

TNTL=T5-T4 57lo 

WRITE ( 6 » 1 6 0 4 ) KH0L,TNTL 5720 

1604 FORMAT (IX, 1 THTIME FqR KInDING »Ho»bh SETS OE/ 5730 

130H NONTOUCrilNG LOOPS* IN SrCONDS, F15.3/) 5740 

5750 
5/6 0 

PROGRAM MA1N--8 5770 

DECODE COMPOSITE SYMBOL CODE b?80 

AND ISOLATE SYMBOLS FROM 5790 

INVERst SYMBOLS 5800 

5Bio 
5820 

NANU = LlL”l 5830 

00 6vl Jl = 1 , NEXP5 5840 

DO 691 J2=l,NTO 5850 

691 POLYU ( Jl , j2) =0 5860 

00 693 Jlsl,NTO ' 5870 
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00 693 J2al,NSPTU 

SFMP0N(J1,J?)=S1aH(D 

SEMP00(JltJ?)=STAR(l) 

SIMB0N(J1,J?)=SH 
693 SIMB00 ( J1 * JZ) =SH 
00 951 J4sl,NT0 
NA(J4)=0 
951 NR(J4)=0 
C 

C DECODE KSORT ( J7 ) AND RECORD TEoMS 
C CONTAINING FEE^DACK SYMBOL *F3* 

JZU=L I L- 1 

00 646 JZ=1,JZU 
KODY=KSOR I ( jZ ) 

1 top ( jz ) =0 

IF (KODY) 715,646,715 

715 CALL DECODE' (KOO t KOD Y* IZ»F B» JZ,SEiMROL* KOOK »KOoI , ITOP, KB ARTS) 
C , 

C ISOLATE NUM. SYMBOLS AND INVERSE SYMBOLS 
C OF KSORT(JZ) 

637 NAK=0 

NAT = 0 ^ 

IF(TZ)646*646»647 
647 CONTINUE 

DO 6^5 NZ=1 , 1Z 
KOZY=KODI (NZ) 

IARGsKODF (N7) 

I F ( I ARG-NRS > 1 340 » 1 34 0 » 1 34 1 

1341 WHITE (6, 13*2) 

1342 FORMAT (] X,3 0 HINCREASE' T He OIMfNSION OF S f AR ) 

1340 CONTINUE 

I F ( KONS ( KOZ y ) ) 65 7 » 657 » 65y 
657 NAK=NAK*l 

IF (NAK-NSPTn-1) 1376, 1375, 1375 

1375 WRITE (6, 137; ) 

1 377 FORMAT ( IX, 4nMNSRl EXCEEDS ' L I M I T- T NCRF. ASE DIMENSIONS , 

1 15HC0NT A INING NSPT) 

1376 CONTINUE 

SIMB0N( JZ*NftK) =5EM30 L (KOzY) 

SEMPON(JZ,N a K)=srAR(IARG) 

NA ( JZ) =NA (J 7 ) *1 
GO TO 645 
659 NATsNAT* 1 

IF (NAT-NSPTu-1) 1381,1380,1380 

1380 WRITF (6, 138?) 

1382 FORMAT ( IX, 4oHnSRT EXCEEDS L I M I T - I NCRF ASE DIMENSIONS , 

] 15HC0NTAINING N$P T ) 

1381 CONTINUE 

SIMBOO ( JZ,NaT) =SEM90L (KUZY) 

SEMPOOt JZ»N aT) =S I AR ( IARG) 

NB ( JZ) =NB ( JZ) ♦ 1 

645 CONTINUE 

646 CONTINUE 


PROGRAM MAIN--9 

St^ARATE POLY InTU ARRAYS FOR THE 
NUMERATOR and DENOMINATOR UF TmE 
TRANSFER FUNCTION 


Sfl 


601 n 
©02f 

603 fUi 
6040 
oOSfpa 
bObrij 
bO 70 
60B0 
609 O 
6 lOti 
6l 1 0 

6 1 2 ft *} 

tnP 
6160 
61.6 

6 1 7 
bl 60 
6l4®~j 

si;iu 

b22o 

6250 
b 2 b 
b27 1 ! 

62 b r 5 

b2 4 o 
630fv 

63 1 
b3? o 
b3 30™ 
6 3 4 1 1 

bisP 

b36 0 
637F| 

6340 

6400. 


0 


641 

642 vJ 

b43n 

644 r | 

b 45 lj 

6460 
04 7 0 

•*"Q 
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c the constant Coefficients in the transfer function 649o 

C ARE SEPARATED INTO ARRAYS FuR THE NUMERATOR 6500 

C AND OENOmInAIOr 6510 

KlKUeKlK-I . ' - — 6520 

00 755 JA=I,MKU b530 

JIB=0 6540 

Jl)=0 6550 

DO 755 JCslfNAWU 6560 

IF (ITOP (JC) ) 753*753,751 6570 

751 JlHsJlB-fl 6580 

POLYU( JA* JlR) =RULY ( JAt JC) 6590 

60 TO 755 6600 

753 JO* Jl) + 1 0610 

POLY (JA,J 0 ) -POLY (JA,JC) ' o6?0 

755 CONTINUE 6630 

CALL SEC0ND'(T6) 6640 

TDF.C0() = T6~Tc; 6650 

WRITE ( 6 » 1 6 0 5 ) 1 DECOD 6660 

1606 FORMAT (1X*36HT1ME FOR DECODING SYMBOLS IN StCONOS , F 1 5« 3/ ) 6670 

6680 
6690 

, PROGRAM MA I N-- 1 0 ' 6700 

make powers of s In output 67m 

transfer function positivl 672o 

" \ o7 30 

b 7 AO 

.MAXIM=0 b75n 

KARU = KlK-l 0760 

DO 522 KAR=i»KAKU 0770 

IF (MS0HT (KAR) )521»522,522 6?80 

521 IF (MAXIM *mSqRT (KAR) )523 « 622 » 5?2 076=0 

523 MAXIm=-msO«T (KAK) b8on 

522 CONTINUE 1 ' 6810 

DO 524 KIT=]»KAHU 6820 

524 MSOR r (KI T) = mAXIM + MS0RT (K l n b8 3o 

6840 

- 6850 

PROGRAM MAlN--.il b860 

PRINT OUT fsjUMERATOR OF b8?o 

The Transfer function ohbo 

0040 
6900 

LUK = 0 . - - 6910 

IKUsLIL- 1 b9?o 

WRITE (6*93!) 6930 

WRITE (6»93(j) bV40 

WRITE <6,92 ^) 6950 

920 FORMAT ( 25^^20 HNUMERA TOR POLYNOMIAL///) 6960 

WRITE <6,92!) 6970 

921 FORMAT ( 1X,0HCDLUMN, 12X,23HSYMbOL FOR GIVEN COLUMN) 6980 

00 905 IKsl.IKU 6990 

IF (ITOP( IK) s 905»905,901 7000 

901 ILUsNA(IK) 70lo 

IF(TLU)710»710»711 7020 

710 ILUsl ' 70 30 

711 JLU*N8(IK) 7040 

IF < JLU) 712»7lH»713 7050 

712 JLUsl 7060 

713 CONTINUE 7070 

LUK=LUK+,i 7080 

WRITE (6,903) LUK, <SIM60n<IK,IL) ♦ SEMPON ( IK, Tl) , 7090 
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1 IL = 1 ♦ ILU) » D /, S H » (S1M8O0UK, JL) , SEMPOD (I K » JL ) ,JL=l,JLU) 

903 FORMAT ( IX, I5*20X,30AJ> 

905 CONTINUE 

WRITE (6.930) 

930 FORMAT (//) 

WRITE (6, 1h2t ) 

1821 FORmAT ( 1 X « 7H POWER) 

WRITE (6,92?) 

922 FORMAT (lX.ttH uF S ♦ 1 7x ♦ 33HCONST ANT COEFS. In THE POLYNOMIAL) 
LML=1 
LMUs7 

IF ( jIB-LMU) 820.618.818 
820 LMIJsJIb 

818 WRITE (6,806) < LO » L0=LML , LMIJ) 

806 FORMAT (2X.7 (OX.bnCOLUMN, I 2) ) 

KROWU=KIk-1 

DO 808 KROWsl .KHOWO 

WRITE ( 6, 81 A ) MSURT (K«OW) , (POLYU (KROW.LM) «LM=LML,LMU) 

810 FORMAT (I5.JH »7<E12.b»4H )) 

808 CONTINUE 

IF (JI8-LMU) 814.814, 812 
812 LML=LML ♦ 7 
LMUsLMU*7 

IF (JIB-LMU) 8l6.8lB.8l8 
816 LMUsJIb 
GO TO 818 
814 CONTINUE 


PROGRAM MA1N--1? 

PRINT OUT DENOMINATOR OF 

The transfer function 


LUK = 0 
IKU=LIL-1 
WRITE (6.93 y) 

931 FORMAT ( 5 OH *<»#**** #■*«##*'»**#**#*•»***»*#«■*«*** ***« #<»<»*#*#*«■#»# ) 
wRITF (6,93n) 

WRITE (6.92J) 

923 FORMAT ( 25X , 22HDEN0M I NA ToR POLYNOMIAL///) 

WRITE (6,924) 

924 FORMAT ( IX.a^COLUMN. 12X.23HSYMB0L FOR GIVEN COLUMN) 

DO 705 I K = 1 , IKU 

I F ( I T OP (IK) , 701.701,705 
701 ILU=NA(IK) 
luk=luk*i 

IF ( ILU) 915.915, 916 

915 1 LU= 1 

916 JLUsNB(lK) 

IF (JLU) 917, 917.916 

917 JLU= 1 

918 CONTINUE 

WRITE (6,703) LUK, (SIMBON ( IK, IL) ?SEMPON( IK, IL) ♦ 

1 1 L= 1 , ILU> , 0 aSH , ( S I M80D ( I ft , JL ) , SEmPOD < I K, JL 5 , JL=1, JLU) 

703 FORMAT ( IX, I5,20X,30A3) 

705 CONTINUE 

WRITE (6,93 0) 
wRITE (6, 1822) 

1822 FORMAT ( IX, 7h POWER) 

WRITE (6,925) 


/bSO 
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925 F0RMAT(1X»8h OF S « 17X» 33HCONSTANT COEF S. IN THE POLYNOMIAL) 7710 

LML=1 7720 

LMU=7 7710 

IF ( JD-LMtJ)S20»518«5l« 7740 

520 LMUbJD 7750 

51B WRITE ( 6 » 506 ) ( LO » LO=LML ♦ LMU ) 7750 

506 FORMAT ( 2X » 7 (8X»6HC0LUMN» 12) ) ' 7770 

KROWU=KIK-l 7780 

00 508 KRUW=1 »KROWU 7790. 

WRITE (6»5lA> MSORT (KROW) , (POl.Y (KROW»LM) ? LmsLML»LMU) 7800 

510 FORMAT ( 1 5 « 3n »7(El2.5»4H )) 7810 

508 CONTINUE 7820 

IF<JD-LMU)5l4 ? t9l*t5l2 - 7830 

512 LMLsLML+7 7840 

LMU=LMU*7 7850 

IF ( JO-LMU) 5l6*518» 518 ' 7860 

516 LMUsJO 7870 

00 TO 518 . 7880 

514 CONTINUE , - 7890 

WRITE (6*930) , . 7900 

CALL SECOND (TEND) \ 7910 

TEXECsTENO-tSTAKT 7920 

WRITE ( 6 » I lb 1 ) TcXEC 7930 

U61 FORMAT (lX«e:7H EXECUTION TIME IN SECONDS » F 1 5 . 3// , 7940 

] H8H AUGU$I i 9 70 VERSION uF SNAP) " 7950 

250 00 TO 1111 7960 

END 7970 

SUBROUTINE sFG(NFIRsT»NLAST.IXPON,WEl0T»SYMBUL»KONSOtMTX»N!tST, 7980 

1LIST»nIn»nOUT ♦ NOU ♦ NQH t L 1 S I S»NDUA,NODB) 7990 

; «0 1 0 
c the following arrays arf associated with the nfjwork b 02 o 

C CHARACTERISTIC NdN (DEFINED IN PROGRAM MAlrW) b030 

DIMENSION JpOw ( 35) t NP(35) »IVV(35) ,NUML(3b) ,ICV(35> ♦1NTPEE(35) 8040 

DIMENSION HNC(35) 8050 

DIMENSION Nf( 35*35) »I8(3b*35) »NS(35»35) 8060 

DIMENSION TYPH(3b) »JO(35) ,L8(3b). ,MSYM(35) b070 

DIMENSION IWUAL(35) »VAL(35) »SYM(35) 80«0 

DIMENSION IoUALX ( 35) »VALX(35> »NUMLX(35) »1NTPE(35) «NOTRFF(^5) 8090 

DIMENSION TYPX(35) *NI)MX(35) ,JBX(35) t L8X(35) ♦SYMX(35) HlOO 

bl20 

C The FOLLOWING ARRAYS are associated WITH THE NFTwORk 8 1 30 

C CHARACTERISTIC N8G bl4o 

DIMENSION Cy AL (100) bl50 

DIMENSION NFlMST (100) *NLAST(100) *IXPON(100 ) »w£I 6T(100) bl6ft 

DIMENSION MARY ( 100) »KONSO( 100) tNEST ( 100) »TYPt( 100) 8170 

DIMENSION SyMHOL ( 100) *MIX ( 100) 8l80 

C««<HH> **##***#*##*#*#*####*#4 tt *«HH» ««• 4MMH» « «■ ««««#«««»«««»««« 8 1 90 

COMMON NF * NS ♦ 18 b200 

COMMON/C2/NNG,NBG , 7 82lo 

' b2?0 

SUBPROGRAM 7A* b230 

DATA Y*G»C» T<J»R»CLtZ/2HY »2H6 »2HC *1H=«2 mr ,2hL «2hZ / 6240 

DATA EtCl»CC»CV*VV»VC/2HE *2HI » 2HCC . 2hC V , ?H VV , 2HVC/ 6250 

DATA Fb/3H Fb/ 1 , b260 

DATA 0NE/3H 1/ 8270 

00 710 IC=1,NNG b2flf> 

DO 710 IK=1»NNG 6290 

NS (ICtlK) =0 , 8300 

710 NF ( I C ♦ I K ) =0 b3l0 
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LINK = 0 

no \52 IG=1,NRG ' 

NEST ( IG) so 
152 KONSO(lG)=0 
WRITE (6,26^) 

260 FORMAT (//) • ' 

wR!TE<6,717) 

717 FORMAT (30X»7RnETwORK/) 

NlAST <1 )*nIn 
1 XPON (15=0 " 

- WFIGT(l)=-l, v .. 

SYMRUL ( 1 ) sFp 

KONSO ( 1 ) =0 

NEST(1)=1 

MOsQ 

LO=0 

L 1 ST= 1 

KLU = 0 

00 5 11=1 »NmG 
INTRF.E ( 1 1 ) *0 
5 JROw ( 1 1 ) =o 
DO 528 1=1 .nOh 

HFAp (5*9) TyPX(I) *NUMX(I)«JHX(1) tLBXd) tSYHX(I)* 
llQUALX(I) «V»LX(1) tNUMLX(l) 

IF(TYPXd) .EQ.COGO TO ljOO 
IF (TYPX( I ) .eQ.CV) 60 TO noo 
IF (TYPX(I) etU.VV) GO ro lJOO 
IF(TYPX(I) .tQ.VOBO TO 1300 
GO TO 1301 

1300 IF { MUML X ( I ) ) 1301* 1302 » 13d 1 
130? WRITE (6, 1303) 

1303 FORMAT ( lX»4yH<H»*ERROR# e *CONTROL SPECIFICATION FOR OFP. SOIJRCF 
17HMISSING) 

NORsO ) 

GO TO 1306 

1301 CONTINUE 
528 CONTINUE 

9 FORMAT (A2*i3,215* 1X»A3* A1,E1?*5*I3) 

WRITE (6*261) 

261 FORMAT (2X*3?HELEMENT ELtwENT INITIAL TtHMlNAL* " 

12RH ELEMENT ELEMENT ELEMtNT NO,) 

WRITE (6* 262} 

262 FORMAT <2X,3qH TYPE NUMbER WORE NODE SYMBOL * 

120H VALUE OF CONTROL) 

DO 601 M=l,NOh 

601 WRITE (6*600)- TTKXtcn*NUMx(M) ♦ JtrXTM) iLBXIM) fSYtfxtMl* - 

1 10UALX (M) » VaLX (M) » NIJMLX ( M ) 

600 FORMAT (4X* A2*6X * 12 ,6X* 1^, 6X» I2*6X* A3* A1 »E12.5,2X * I?) 

CALL FTREE(tYPX* JHX,LBX» inTRE»NOTREF*NOO»NOR) --- - 

SUBPROGRAM tHt 

WRITE (6*5U) ■ ' 

518 FORMAT (30X* 13HTREE SELEClEO) 

NUMU=NOD-l 
00 21 NUs 1 *NUMU 
1 0= I NTRE ( NU ) 

NUMC = NUMX ( 1 0 ) 

TYPB (NUMC)=TYPX ( 10) 

JB (NUMC) =JtiX ( 10) 

L8 (NUMC) sLb x ( TO) 

SYM (NUMC) =SyMX ( 10) 


8320 fl 
8 3 3 a U 
83 AO 
8350 rn 
83601 ( 
8370 U 
8380 
83900 
840 oLi 

b A 1 0 
842 OrT 
84.3 0 | 
b440^ 
8450 
8460p| 
84T0-J 
8480 

849 Op, 

850 C* j 
8510^ 
85?0 

853 cn 

854 qj 
8550 
8560™ 
857 f I 
858fU 
bBRfi 
dbO'H 

86,1 (J 
8620 
b6 3 0 

864 ( J 
b6SiLI 
8660 
86 7 Cl 

o6Q ii 

8690 
8/0 a 
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IOUAL(NUMC) alOUALX (10) 

VAL(NUMC)=V a LX(I0> 

NUML (NUMC) =nUmLX (10) 

INTREE (NUMC) =1 

WRITE ( e> 9 5 1 7 > TYPB (NUMC) « NuMC ♦ JH ( NUMC ) * LB (NUMC) *SYM(NlJMC) » 
1IQUAL (NUMC) ,VAL (NUMC) «NUmL (NUmC) 

517 FORMAT (4X»A2»6X»l2,6X«l2,6X*T2»6X»A3»Al»El2,5 t 2X»l2> 
KIU=KLU*1 
LINC (NUMC) =0 

IF (TYPB (NUM9) .nE. W) GO 10 3 
MO=MO*l 

IVV(mO)=nUMC s 

3 IF (TYPb (NUMC) .NheCV) GO 1 0 A 

lo=lo*i <■ 

ICV(LO)=NUMC N 

4 JF s jh (NUmC) 

LFsLB(NUMC) 

I B ( JF«LF) sNiiMC 

IB(UF« JF) aNiiMC - 

JROW ( JF) ajRoW ( JF) ♦ 1 
JROJ=JHQW ( JT ) 

NF ( JF ? JRO J ) S LF - 

NS(JF ? LF)sl 
JROW (LF) sjRQW (LF) +1 

JROLsjROW (LF ) 

NF ( LF « JROL ) = JF 
NS(LF,JF)=-1 

21 CONTINUE 

WRITE ( 6 ♦ 260 ) 

wRITE(6,715) 

(15 FORMAT (30XOHSFG/) 

DO 13 I L L — 1 » NOO 
JR0T=JRUW ( ILL) >1 
13 NF ( ILL « JROl ) =0 

SUBPROGRAM FCF 

THIS PROGRflM GENERATES SIGNAL FLOW GRAPH INFO. 

FROM BRANCH NOUt TO LINK NODE 
NO0Y=NOB 
151 CONTINUE 

NFS = 0 
LON = 0 

IF (KLU-NUB) 02*360 *532 
532 L INKsL 1 NK* 1 

1F(N0TREE(LINK) )534,534»b32 
534 NUMCsNUMX (LINK) 

TYPE (NUMC) -TYPX (LINK) 

JK=JBX(LINK) 

LKsLH X (LINK) 

SYM(NIJMC) =STMX(LINK) 

1QIJAL (NUMC) =10UaLX (LINK) 

CVAL (NUMC) =VALX (LINK) 

NUM9=NUMLX (LINK) 

TYP2 = TYPE(N||MC) 

CVALUsCVAL (NUMC) 

KLU=KLU«1 

LINC(NUMC)=1 

KOEPSsO 

KANSO=0 

IF ( Type (NUMC) eEO.CDGO TO 117 
1F(TYPE(NUMC) .e«cG)GO TO 119 


6930 
b940 
6950 
6960 
6970 
6980 
8990 
9000 
9010 
9020 
9030 
9040 
9050 
9060 
9070 
9080 
9090 
9100 
9110 
9120 
9130 
9140 
9150 
9160 
9170 
9 180 
9190 
9200 
921 0 
9220 
923/t 
9240 
9250 
926o 
9270 
9280 
9290 
9300 
9310 
9320 
9330 
9340 
93So 
9360 
9370 
9380 
9390 
9400 
9410 
9420 
9*3p 
944 0 
946o 
9460 
9470 
9480 
9490 
9500 
9510 
95?0 
9530 
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1F(TYPE(NUMC) .EO.Y) GO TO 119 
IF (TYPE (NUMC> ,EO.R) GO TO 700 
IF (TYPE(NUMC) «EU.Z) GO TO 700 

if ( Type (numc> .to. a go to l ? i > 

KDEPS=1 

IF (TYPE (NUMc> .EU.E) GO 10 123 
IF (TYPE (NUM r ) ,EU.CI)GO To 123 
IF (Type (numc) .Eu.vo ho T o 165 

IF ( TYPE (NUMr ) ,E(J.CC)GO To 265 
117 I XPS = - 1 
KANSO=l 

GO TO 123 ' 

1 19 IXPSsO 

GO To 123 
700 I XPS=0 
KANS<)=1 
GO TO 123 
121 IXPS=1 

123 CALL TREP(Jk*LK*NF*NH*NPl) 

ifin=numc 

U9 LON=LON+1 
NP1=NP(L0N) 

NP2=NP ( LON+ 1 ) 

107 INI T-IB(NP1 ,NP2) 

109 $IGN=NS (NPl ,NP2> 

IF ( K EPS) I67 * 167*169 
167 IF (IOUAL(NUmC) .EQ.IQIGO TO 111 
NES=1 

CONSTsSIGN 
GO T() 126 

111 const=sign*cvalu 
125 LIST*LIST*1 

IF (NES) 502*503*502 ' 

502 NFST(LIST)=1 

503 KONSO(LIST) =NANSO 
NF IH5T (LIST) =INIT 
NLAST (LIST) =1FIn 
SYMBUL (LIST , = $YM ( IF IN) 

IXPON (LIST) = 1 XPS 

IF ( KONSO ( L I ST ) )505*505*50A 

504 WEIGT ( L 1 ST ) =1 . /CONST 
GO TO 506 

505 WEIGT (LIST ) = CONST 

506 MAPY (NUMC) =i'lST 
127 FORMAT (3I5,El2.5) 

129 FORMAT ( A4 ) 

SUBPROGRAM *0* 

THIS PROGR/\" GENERATES SIGNAL FLOW GRAPH INFO, 
FROM LINK NODE TO BRANCH NOOF 
169 CONTINUE 

IF (TYPb( INI T > .EU.E) GO TO 201 
IF ( T YPB ( INI T ) .EU.CDbO To 201 
IF (TYPB ( InIt) ,E0.VV)G0 To 201 
IF (TYPfcU INIT) .EQ.OMgO To 201 
LTST=LIST<-1 

IF (TYPM ( InIt 1 • E 0 • R ) G 0 To 133 
IF (TYPH ( InIt) ,E0.2>G0 TO 133 
IF (TYPB ( INIT) .EU.G) GO TO 702 
IF (TYPb(lNlT) ,EO.Y)GO TO 70? 

IF ( TYPB ( iNlj) .EN.CL) GO To 135 


y7 7 fU 

9730 



9B6'! 
98 7 <r 
9 HB 0 
9 BwH 


1 0 1 f? 
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IF (TYPB(InIT) .EU.C1G0 TO 137 10150 

133 IXPOn(LIST)=0 10160 

60 TO HI 10170 

702 IXPON (LIST) sO " lOlflO 

KONSO (LIST) si 10190 

GO TO HI 10200 

135 TXPON(LIST) =1 10210 

GO TO 141 10220 

137 IXPON(LISF) =-l 10230 

KONSO (L 1ST ) = 1 10240 

141 IF ( IUUAL ( INI T ) •£()• IQ) GO 10 139 10250 

NEST(LlST)=i 10260 

WFIGT (LIST) _»1 .»SIGN 10?7o 

GO TO 147 . ' 10280 

139 IF (KONSO (LUl )) 608,608, 607 10290 

607 wFTGT (LIST) ="SIGN/VAL ( INIT) ' 10300 

GO TO 147 10310 

608 WFIGT (LIST) =-SIGN»VAL('INlT) 10320 

147 NFIWST (LIST ) = I F 1 N 10330 

NLAST (USD »INI T - 10340 

SVH8UL ( LIST ) =STM ( I NT T) 10380 

201 NPLA = NPL-l“i-ON 10360 

IF (NPLA) lbl , 151 ♦ 149 1037O 

10380 

SUBPROGRAM E 10390 

this program sets up sfg info, for vc 10400 

Type C0N1H0L SOURCES 10410 

165 NUNO=NUMB 1042n 

IF(TNTHEE(NuMH) ) 163,163.161 10430 

163 LIST=LIST*1 10440 

NFIRST (LIST ) *NUmh 10450 

N0BY=N0bY*1 10460 

NLAST (LIST l-NOHY 10470 

SYMBUL (LIST) =SYM (NiJMH) 1 10480 

NUN0=NObY 10490 

IF (TYPE (NUMw) .EQ. Y) GO TO 912 10500 

IF ( TYPE (NUMq) # E(J.G) GO TO 912 ' 10S)o 

IF (TYPE (NUMa) .EO.OGO TO 914 10520 

IF (TYPE (NUMq) .EU.CLIGO TO 916 n 10530 

NUNO=0 10540 

IXPON (LIST ) =0 - * 10550 

GO TO 918 . 1056p 

912 IXPON (LIST) =0 10570 

KUN0=1 10580 

GO TO 918 10590 

914 IXPON(LIST) =-l * 10600 

KUN0=1 1 061 0 

GO TO 918 10620 

916 IXPON (LIST ) :1 ' 10630 

KUNO = 0 10640 

918 IF ( I QUAL (NUmU 1 • 1 0 • I Q ) GO TO 92(1 10650 

NEST (LIST) — 1 1066 n 

WEIGT (LIST) 41. IO 670 

GO TO 209 lOSHo 

920 IF ( KUNO ) 922 » 922 t 924 10690 

922 WEIGT (LIST) -CVAL (NUMB) 10700 

GO TO 209 ] 0 7 1 0 

924 hE IGT ( L 1ST )= 1 . /CVAL (NUMB) 107?o 

209 KONSO (LIST ) =1 10730 

161 L I ST = L 1 ST + 1 10740 

NFIRST (LIST) =NUNO lu780 
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NLAST (LIST) =NUMC 
SYMRUL(LIST) aSYM(NJMC) 

UPON (LIST) =0 

IF(IQUAL(NUmC) .t(J.IO)GO TO 1 7 1 
NEST (LlST)=i 
WEIGT (LIST) -1. 

GO TO 203 

171 WEIGT (LIST) aCVALU 
203 CONTINUE 
GO TO 123 

SUBPROGRAM 

THIS PROOR/lM SETS UP SFu INFO. FOR CC 
TYPE CONTROL sources 
265 HUN0=NUMB 

IF ( TNTHEE (NiiMB) ) 621 ,t>21 .620 

620 LIST=LIST*1 
NFIRST (LIST ) sNUMB 
NOBY=NOBY«1 
NLAST (LIST ) sNOBY 
SYMBOL ( L ISl )=SYM(NJMB) 

MUNO=NOBY 

IF (TYPB (NUMr) .EU eZ) GO 1 0 233 
IF (TYPH(NUMw) .EU.R) GO TO 233 
IF (TYPB (MUMr) .EO.CDGO To 235 
I F ( T YPB ( NUMo ) • E<J • C ) GO TO 237 
KUNQ=0 

IXP0N(LIS( ) aO 
GO TO 24). 

233 IXP0N(LIS1 )=U 
KUNO=l 
GO TO 241 

235 IXPON (LIST) =-l 
KUN0=1 
GO TO 241 
237 IXP0N(L1ST ) =1 
KUNO=0 

241 IF ( IOUAL (NUvtb) ,EQ. Io) GO fO 230 
NE$T(LlST)=l 
wEIGT (LIST ) =1 
GO TO 247 

239 IF (KUNO) 900,900, 902 
900 WEIGT (LIST) rVAL(NUMR) 

GO TO 247 

902 wEIGT (LIST ) =1 ./VAL (NUMB) 

247 KONSO ( L I ST ) - 1 • 

621 LIST=LlST*l 
NFIRST (LIST) =MUNO 
NLAST(LIST)=NUMC 
SYMBUL (LIST ) =SYM (NJMC) 

UPON (LIST) -0 

IF(TQUal(nUmC) . to. IQ) GO TO 27]' 

NFST (LIST) — 1 
WEIGT (LIST).. 1. 

GO TO 281 

271 WEIGT (LISD-CVALU 
GO TO 281 
281 CONTINUE 
GO TO 123 


SUBPROGRAM #0* 
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N0BY=N06Y* 1 11980 

NLAST (LISD =NUBY 11990 

SYMBUL (LIST) ®SYM (NUNO) 120QO 

IF (TYPB (NUNo) .EO.Z) GO TO 433 12010 

IF (TYPB (NUN n ) .tU.H) GO TO 433 12020 

IF (TYPB(NUNo> .EUeCLJbO To 435 12030 

IF (TYPB (NUNn) .EU. C) GO 10 437 12040 

K U N 0 s 0 12050 

1XP0N (LIST) =0 12050 

GO TO 441 12070 

433 IXPOM(LIST) =U 12080 

KUN0=1 12090 

GO TO 441 12100 

435 IXPON(LIST) =1 12110 

KUNO=l 12120 

GO TO 441 121.3" 

437 IXPoN(LIST) =-l > 12140 

KUNOsO 12150 

441 IF ( tQUAL (nU n O) .tU.KDGO 10 439 ... 12160 

NEST(LIST)=1 ■ ' ' 12170 

WFIST (LIST ) si, ' 1 2 1 HO 

GO TO 448 1219" 

439 IF ( kUNO) 902 *908*910 12200 

908 WFIGT (LIST) =VAL (NUNO) 12210 

GO TO 448 12220 

910 WEIGT (LIST) =l./VAL(NUNO) 12230 

448 KONSO (L 1ST) =1 12240 

447 CONTINUE - 12250 

NUN0=N08Y 12260 

461 L I ST=L 1ST + 1 12270 

NFIRST (LIST) =NUNO / 12280 

NLAST (LIST ) =Lt 12290 

SYM8UL (LIST ) =SYM(LI) 12300 

IXPON(LISf) =0 12310 

IFdQUAL(Ll) *FO.IQ)GO TO 471 12320 

NEST ( L I S T ) s j 12330 

WEIGT (Lisr)-l. . 12340 

GO TO 403 12350 

471 WEIGT (LIST) sVAL(Ll) 12360 

403 CONTINUE 12370 

405 CONTINUE 12380 

12390 

SUBPROGRAM ' 12400 

GENERATING OUTPUT NODE OF SFg 12410 

12420 

515 CONTINUE ' • 12430 

IF (NOUT) 514,512*514 12440 

612 CALL THEP(N n UA,NODB»NF*NP*NPL) 12450 

NOUTsNOBY* 1 ■ , r 12460 

MOPU=NPL~ 1 12470 

00 510 MOpsftMoPU 12480 

NlsNP(MOP) 12490 

N2«NP (MOP+l) 12500 

LIST=LIST*1 125lo 

NFIRST (LIST) =IB(N1 t N2) ' 12520 

NLAST (LIST) sNOUT 12530 

SYMBOL (LIST) =ONt 12540 

I XPON (LIST) = 0 12550 

KONSO (LIST ) =0 12560 

NF.ST ( L I ST > =0 12570 

510 wEIGT (LIST )=NS(N1*N2) 1 - 125mo 
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C THIS PROGRAM SETS UP SFg Info, for VO 

C TYPE CONTROL SOURCES 

160 IF ! MO) 460 »4fc>0 . 3b4 
364 1)0 305 MI = 1 .mo 
KI*I VV (MI ) 

NUNOsNUML (Ky ) 

IF (LINC (NUNU) ) 361 » 361 » 363 
363 LIST*LIST.1 

NF1RST ( L I S 1 ) »NUML(KI) 

N0BY=N08 Y+ 1 
NLAST (LIST) sNUHY 
SYM8UL (LIST ) *SYM (NUNO) 

IF (Type (NUN n) ,KU. Y) GO TO 333 
IF (TYPE (NUN n > .EO.G) GO 10 333 
IF(TYPE(NUNo> .E8.OG0 TO 335 
IF(TYPElNUNn) .EU.CUGO To 337 
KIJNO = 0 

IXPON (LIST) sO 
GO TO 341 

333 IXPON (LIST) =0 
KUN0=1 
GO TO 341 

335 IXPON(LIST) =**1 
KUN0=1 
GO TO 341 
337 IXPON (L 1ST ) =1 
KUNOsO 

341 IF ( IQUAL (NUnO) *EW. IQ) GO TO 339 
NEST (LIST) =1 
WFI6T (LIST) =1. 

GO TO 34ft 

339 IF (KIJNO) 904,904.906 
904 WFIoT (LIST) =CVAL (NUNO) 

GO TO 348 

906 WEIGT (LIST) al ./CVAL (NUNO) 

348 KONSO ( L I S T ) - 1 
347 CONTINUE 
NUNOsNOBY 
361 LIST*LlST*l 

NFIRST (LIST) =NUNO 
NLAST (LIST) -.KI 
SYMBUL(LIST)=SYM(KI) 

IXPON(LIST) =0 

IFdOUAL(Kl) *EU.IQ)GO TO 371 
NEST (LIST)*] 
wFIGf (LIST ) -I. 

GO TO 303 

371 WFIGT (LIST) *VAL(K1) 

303 CONTINUE 
305 CONTINUE 

SUBPROGRAM *rl5* 

THIS PROGRAM SETS UP SFg InFO. FOR CV 
TYPE CONTROL SOURCES 
460 IF (LO) 515.5ib.464 
464 DO 405 MI=l,LO 
Ll*ICV (MI ) 

NUN0=NUML (LI ) 

I F ( L I NC ( NUN() ) ) 46 3 * 463 « 4b l 
463 LIST = LIST<*1 

NFIRST(LIST)=NUML(LI) 
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511 CONTINUE -- - 12590 

514 NFIRST ( 1 ) sN 0 Ur 12600 

48? IF (LlST6)48ft»486»1200 12610 

1?00 wRITfc(b,263) 1 2620 

263 FORMAT (1X,37HINITIAL TERMINAL EXPONENT BRANCH » 12630 

135HRNANCH i IF SYMBOL 1 If- SYMBOL) 12640 

WRITE <t>»264) - - - 12650 

264 F ORMAT ( 1 X * 3 JH NODE NODE OF S VALUE.* 12660 

137H SYMBOL IS iNVERTtD IS USED) ' 12670 

DO 1202 Js l , L I S I 12680 

WRITE (6,46c;) NFIRST(J) *NLAST(J) ,IXP0N(J) fWFIGT (.J) , 12690 

1SYMBUL ( J) »KONSO( J) ,NES1 ( J) 12700 

4H5 F0RMAT(3X»I?»7X,I2,6X*I2»4X»E12.15,1X,A3,8X,I2,14X,I2) 12710 

1202 CONTINUE % 12720 

466 CONTINUE 12730 

v 12740 

SUBPROGRAM FJ * ■ 1.2750 

This PROGRAM ORDERS SPG INFORMATION 12760 

FOR INPU1 TO MAIN PROGRAM 12770 

DO 67 J=1 ,NoG ' 12780 

87 MIX(J)=J « 12790 

KONUsLIST-1 12800 

DO AO KON=l,KONU 12810 

ILJsKONw 1 , 12820 

ILsKON 12830 

GO TO 83 12840 

81 MXL = m I X ( I L ) 12850 

MIX ( IL) =MIX ( IU) , , 1286a 

MIX ( IU) =MXL . 12870 

ILsIL-l 12880 

llJsIU-1 12890 

IF(IL) 80,80,83 12900 

83 ■ MIUSMIA(IU). \ 12910 

MIL s M I X ( IL) . ' “ 12920 

IF (NFIRST (MlU) -NFIRST (MIL) ) 81 , 8y, 80 12930 

84 MXL S M I X ( IL ) 12940 

MIX ( IL) =MIX ( IU) 12950 

MIXUU)=MXL 12960 

IL-IL-1 12970 

IU=IU-1 12980 

IF (IL) 80,80,8? • , 12990 

82 MIUsMlX(lU) 13000 

MILsMIX < IL) - ■ ’ ■ 1301n 

IF (NFIRST (MiU)-NFIRST (MIL) ) 80,b9,80 13020 

89 IF(NLAST (MlU) -NLAST (MIL) ) 80*80*04 1J03O 

80 CONTINUE 13040 

1305 CONTINUE 13050 

RETURN . 1306O 

END . 13070 

SUBROUTINE fTR£E(TYPX»JBa,L3X,1NTRE,NOIRP. E,NUU,N0H) 13080 

13100 

C The FOLLOWING ARRAYS ARE ASSOCIATED WITH THE NETWORK 13110 

C CHARACTERISTICS NUN* AND NSpT (OEMNEU IN PHOgPaM MAIN-1) 13120 

DIMENSION TyPX(35) ,JBX(35) *LHX(35) ,INTRE(35) *NOTPEE(35) 13130 

DIMENSION N r ; ( 35) * NF (39*35) *KCOL (20) 13140 

CtttttHHHHHttttftfttHUHl****************************************************** 13150 

COMMON/ C2/NnU,NBG ' 13160 

DATA E 9 V V , C V /2HE *2 HVV*2hCV/ 13170 

DATA R*CL*C,Y, Z/2HR ,2ML *2HC »2HY ,2HZ / )31ro 

DATA 6/2HG / 13190 
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DO 40 12=1 »NNG 
DO AO 13*1 *nNG 
40 NF ( 1 2 » 1 3 ) =0 
M*0 
K = 0 
KC = 0 

DO 1 T=l»NOn 
1 KCOL ( I ) = 0 
DO 3 T7=1 ,NqB 

3 NOTRFF(I7)=c 
1=0 

5 1*1*1 

I F ( T YPX ( I ) .eO.EIGO TO 10 

6 I F ( T YPX ( 1 ) ejrO.VVlGO TO 1() 

8 IF ( TYPX ( I) «eO.CV)G 0 TO lu 

00 TO 4 
10 K=K*1 
14 1NTRE (K) =1 
JWXlsJBX ( I ) 

KCOL ( JBX1 ) =kCOI. < JBX 1 ) +1 
KCOL 1 =KC0L 1 |bX 1 > 
mF(J8X1«KC0lD=LBX(I) 

IBX1=LBX ( 1 ) 

KCOL(I«X1)=kCOL(IBX1)*1 

kcol2=k:col ( rbx i ) 

NiF (IBXltKCOi-2)=JKXl 

NOTREE U ) =1 
1F(K-N0D«1) ? »22»22 
? I F ( M ) 4 » 4 9 1 2 

4 IF ( I “NOB ) $ * 1 2 » 1 2 - 

1? m = M «■ 1 

IF (TYPX (M) .EO.H) 00 TO lb ' . 

IF (TYPX (M) .pU.G) 00 TO lb 

17 IF (TYPX (M) .eU.CLIOO TO lb 

18 IF (TYPX (M) opU.C) 60 TO 16 
10 IF (TYPX(M) « 1 0 , Y ) 6 0 TO lb 

20 IF (TYPX ( M ) .E'O.Z) GO TO lb 
IF (M-NOB) 12,22^22 

16 NINX=JBX (M) 

NODI X=LBX (M) 

CALL TREP(N1NX*N0UTX»NF ♦nP,NPl> 

I F ( NPL ) 2l*2i»12 

21 T=M 

60 TO 10 

22 CONTINUE 
return 

END 

subrout I me tbep (NIM»N0UT »NF*NP»NPL) 

c ft ft ft ftftft ft* * ft ft * tt ft ft ftftft ftftftftftftfto ft ft ttftftftftftftftftftftftftftft ftftft ftftft ft ftftft ft ft ft ftft ftftft ft ft «» ft 


C The FOLLOWING ARRAYS are ASSOCIATED with The NF1W0RK 
C CHARACTERISTIC ( oEF I NED IN PROORAM MAIN-1) 

DIMENSION Jx (35) ,NP(3b) » JMEM(3b) ,KMFM<35) 

DIMENSION NF ( 3 S * 3 5 ) 

c ftftft ft »**ft*ft»*ftftft« ftftftftftft ft ftftft ftftft* ft ftftftftftftftftftftftftftftft ft ftftft** «««« ft ft* ftftft ft ft 9 ft ft* ft ft 

C0MM0N/C2/N n G*NBG 
DO 80 15=1 »mNG 

JX ( 15) =0 - 

NP ( 15) =0 
JM£M 1 5 ) =0 

80 KMEM(I5)=0 - 
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NPL = 0 13810, 

JX ( 1 ) sN I N 13B20 

JX ( 2 ) sN I N 13830 

1=1 13840 

J=NIN 13050 

NP ( I ) =N I N 13860 

20 K=0 13870 

25 K = K+1 13880 

IK(MF( J*K)-mOUT)30?50*30 13890 

30 IF (NF(J,K) >34,32,34 13900 

32 IF ( J-NlN) 60, 100,60 13910 

13920 

FLOWER CHECK 13930 

34 NJKsmE(J,K) 13940 

IF (NJK-JX ( I )) 45,25,45 13950 

13960 

STORE AND REMEMBER VERTEX 13970 

45 1=1+1 13980 

NP(I)=NF (J»K> 13990 

JMEM(I)=J 14000 

I A= I + 1 14010 

JX ( I A > =NF (J,K) 1 40?o 

42 J=NF ( J , K ) . \ 14030 

KMFM(1)=K 1 4040 

GO TO 20 1 4 050 

14060 

BACKSTEP 14070 

60 JaJMEM(I) 14080 

K=KMEM(I) 14090 

1=1“ 1 14100 

GO TO 25 - 1411 r, 

14120 

FINAL PATH VERTEX AND PaTH LENGTH ' 14130 

50 11 = 1 + 1 - - - 1 *+ 1 40 

NP(II)*NOUT ' 14150 

62 NPL=II , 1 4 1 6 0 

100 CONTINUE 7 - - 14170“ 

RETURN 1 4 1 8 0 

END \ I4l9n 

SUBROUTINE aRRAY ( JSIG,XC0N, JXPO, JKOD, POLY, LIL,KIK) 14200 

C#####*##**##**#####** ###4*#*#*#**#***#####*##*#*## **##****#*#**»***#*# 142 10 

14220 

C THE FOLLOWING ARRAYS ARE ASSOCIATED WITH THE NF I WORK 14230 

C CHARACTERISTICS NTO, AND NExPS (DtFlNE'D IN PROGRAM MAIN-1) 14240 

DIMENSION MSUWT ( 15) ,KSORT (150) »P0LY ( 15, ISO) 14250 

14260 

COMMON/C1/MSORT ,KS0Rl 14270 

’ COMmON/C3/Ne6PS*NTO 14280 

MMX=0 ; 1 4290 

NNX=0 14300 

IE (KlK-1 ) 3»J?2 , 3 14310 

3 MMUsK I K- 1 1432n 

DO 2 MM=1,MmU ' 14330 

MMX = MMX + 1 14340 

IF ( JXPO-MSOrI (MM) ) 2, 10,2 ’ v 14350 

2 CONTINUE 14360 

22 MSORT (K1k)=jXPO 14370 

MMX=KIK 143ro 

K I K = K 1 K + 1 14390 

IF(KIK-NEXPs-l) 1386,1385,1385 14400 

1385 WRITE (6, 1307) ' 144ln 



-123- 


1387 FORMAT ( IX* 42HS-F0WETR EXCEEDS L 1 MI T- J INCREASE DIMENSIONS * 

116HC0NT AlNiNG NEXPS) 

1386 CONTINUE 

10 IK (L IL“1 ) 11*24,11 

11 NNUsLIL -1 

00 12 NN = l ♦ imNU 
NNX=NNX*l 

IF (JKOU~KSOrT (NN) ) 12,20, 12 

12 CONTINUE 

24 KSORT ( L I L ) = JKOO 
NNX = L IL 
LIL=LIL+1 

IF (LlL-NTO-i) 1367*1365,1365 

1365 WRITE (6, 1366) 

1366 FORMAT ( IX, 4*HnO. OF TERMS IN OUTPUT F.XCEtOS L 1 Ml T - I NCRF ASP » 
125HDIMENSI0NS CONTAINING NlO) 

1 367 CONTINUE 

20 POLY (MMX, NNX) =PULV (MMX, NnX) - fXCON* (-1 ,) **JSlfi 
RETURN 
ENO 

SUhROUT INE Ot CODE (K0O*K0uY» 1 2 , F8 , J2 , SEMBOL , KuDF , KOU I , I TOP , KR AS I S ) 
C^o***^##******^#********^'!***** *********** a******************,,********** 


C The following ARRAYS ARE ASSOCIATED WITH The network 
C CHARACTERISTICS NSP f » AND NlO (DEFINED IN PROGRAM MAlN-1) 

DIMENSION I T OP ( ISO) , SEMBol. (20) *KODF (20) ,l\0Dl (20) 

C****«*«*<H»******.&4<HKHMH,fl.»0<H><»**«>**«tt**<Hy******ft<H»**a,HXH>*<><J«*<HHHM»*** 

C0Mm0N/C4/NsRT 

17=0 

M=K8ASIS“I 
DO 3 J=l,KOn 
IPOWEP=M. ANU.KOUY 
IF ( I POWER ) 3,3,2 

2 IF ( SEMBOL ( J ) • EO, F 8 ) 60 TO 4 
17=12*1 

IF (12-NSPT-i ) 1371 , 1370, 1 370 

1370 WRITE (6, 1372) 

1372 FORMAT ( IX, 4 o Hn 0« OF SYMBOLS PER TERM EXCtEDS OUTPUT- I NCRF « SE , 
126H DIMENSIONS CONTAINING NSPT) 

1371 CONTINUE 

KODF ( 12) sIPqWER 
KODI (IZ)=J 
GO TO 3 
4 I TOP ( J2 ) = 1 

3 KODY=KOOY/KrASIS 
RFTURN 

END 
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